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INTRODUCTION 


I 

A.  IMPORTANCE  OF  DETERMINING  OCEAN  THERMAL  STRUCTURE 

An  important  physical  property  and  product  cf  various 
|  dynamical  processes  occurring  in  the  oceans  is  the  thermal 

structure.  Although  it  is  derived  from  only  a  few  mecha¬ 
nisms  such  as  air-sea  thermal  energy  exchange,  forced  ar.d 
I  free  convective  mixing  and  horizontal  ocean  thermal  advec- 

tion,  the  interaction  of  these  processes  along  with  ether 
hydrodynamic  motions  can  become  immensely  complicated. 
Since  the  temperature  profile  has  important  implications  for 
climatic  weather  prediction,  commercial  fisheries  and  mili¬ 
tary  acoustic  surveillance,  to  mention  a  few,  the  ability 
to  predict  this  structure  is  of  major  importance. 

Climatological  atlases  present  an  averaged  profile  for  a 
certain  period  and  location.  However,  this  results  in  a 
generalized  thermal  structure  which  is  not  accurate  enough 
to  accomplish  the  precise  needs  cf  underwater  sound  propaga¬ 
tion  forecasting.  Recent  developments  in  numerical  modeling 
of  the  upper-ocean  temperatures  have  produced  vast  improve¬ 
ments  in  the  analysis  and  prediction  of  mixed  layers  and 
other  synoptic  oceanic  features  (Clancy  et  al. ,  1981).  It 


is  inevitable  that  this  method  will  reach  the  sophistication 


cf  today's  meteorological  forecasting,  but  with  the  complex¬ 
ity  of  the  physical  processes  involved  and  with  present  data 
sparsity  ir  many  ocean  areas,  many  inaccuracies  still  exist. 

The  bathythermograph  (BT)  provides  an  accurate 
temperature  profile  at  a  point  location  in  the  upper  ocean. 
Since  the  vertical  thermal  structure  over  a  large  area  pro¬ 
vides  the  main  factor  which  affects  the  propagation  charac¬ 
teristics  cf  underwater  acoustics,  the  need  for  many 
temperature  profiles  in  Havy  ASH  operations  is  obvious.  The 
expense  involved  in  ship  and  aircraft  deployment  of  BT's  for 
an  ocean  area  of  interest  often  prohibits  the  use  of  this 
most  accurate  means  of  obtaining  a  detailed  thermal 
descripticr. 

B.  SATELLITE  DETERMINATION  OF  SUBSURFACE  THERMAL  STRUCTURE 

The  study  of  the  oceans  by  satellites  aas  become  a  major 
arena  for  scientific  scrutiny  and  investigation.  Already  it 
has  given  added  Insight  into  the  numerous  processes  that 
govern  oceanic  behavior.  The  Gulf  Stream  eddies  and  fronts, 
the  coastal  upwelling  off  the  Somali  and  Arabian  coasts,  and 
other  applications  in  biological,  geophysical  and  dynamical 
fields  are  becoming  well  Jcncwn.  However,  in  spits  of  great 
progress,  satellite  remote  sensing  techniques  for  observing 
oceanic  processes  are  far  from  being  fully  realized. 
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The  fundamental  problem  which  continuously  arises, 
besides  atmospheric  effects  which  will  not  be  discussed,  is 
that  satellite  observations  are  made  at  the  surface  alone 
and  that  all  conclusions  about  subsurface  processes  are 
obtained  from  inference  rather  than  by  direct  measurement. 
Consequently,  the  applications  of  satellites  for  oceanogra¬ 
phy  are  limited  to  those  phenomena  that  generate  surface 
signatures  detectable  by  electromagnetic  waves.  Once  these 
surface  manifestations  are  observed,  it  is  necessary  to  have 
knowledge  of  the  physical  processes  involved  to  interpret 
their  meaning.  Currently  there  are  few  methods  available 
where  a  qualitative  "feel"  for  subsurface  features  can  be 
directly  expressed  as  empirical  relationships  from  remotely 
sensed  measurements  at  the  sea  surface.  If  a  linkage  is 
found  to  exist  between  the  air-sea  interface  and  the  thermal 
structure  below,  a  satellite  would  then  be  the  perfect  plat¬ 
form  to  collect  sea  surface  data  rapidly  over  all  the 
world's  oceans.  This  surface  information  could  then  be  used 
to  infer  the  subsurface  vertical  temperature  structure. 

C.  EOHPOSE  AEE  GCALS  OF  TEESIS 

The  primary  purpose  fcr  the  ongoing  research  into  the 
or  edict  ion  of  subsurface  thermal  structure  is  to  understand 


the  processes  and  relationships  that  govern  the  ocean's 
behavior.  Holes  that  are  performed  by  ships  and  aircraft 
(via  expendable  bathythermographs)  in  determining  the 
temperature  profiles  could  be  played  better  by  satellites, 
provided  a  "link"  is  found  to  read  the  below-surf  ace  events 
from  the  sea  surface  signatures.  Therefore,  if  strong  cor¬ 
relations  and  relationships  could  be  ascertained  between 
satellite-derived  sea  surface  temperatures  and  vertical 
temperature  profiles,  then  the  areas  of  acoustical  oceanog¬ 
raphy  and  naval  tactical  applications  would  benefit  greatly. 

The  goal  of  this  thesis  is  to  determine  the  possible 
correlations  that  exist  between  the  subsurface  thermal 
structure  and  the  temperature  at  the  surface.  Then,  from 
these  relationships  and  use  cf  regression  formulas, 
predicted  vertical  temperature  profiles  could  be  computed 
from  an  input  cf  sea  surface  temperature  alone. 
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II. 


CHARACTERISTICS  OF 

pscirurrreCT - 


A.  BIND  SISTESS  AND  CURRENTS 

The  Northeast  Pacific  has  been  studied  extensively  for 
many  years  and  it  has  revealed  a  complex  arrangement  of 
environmental  conditions.  Severe  storms,  high  winds,  varied 
water  masses,  numerous  currents  and  unique  temperature  and 
salinity  structures  are  a  normal  occurrence  for  the  region. 

In  the  winter  the  Ncrth  Pacific  is  usually  under  the 
influence  of  the  Aleutian  Low  in  the  atmosphere  while  for 
the  summer  mentis  a  high  pressure  system  dominates  (Uda, 
1963) .  The  distribution  of  the  mean  surface  winds  are  con¬ 
trolled  by  these  two  systems.  At  the  location  cf  the  new 
retired  Ocean  Heather  Station  ”PW  (OHS-P,  positioned  at  SON, 
145H  in  the  Gulf  of  Alaska)  the  prevailing  winds  are  from 
the  southwest  (Tafcata,  1S65)  .  From  the  influence  of  these 
wind  patterns,  the  surface  waters  flow  eastward  across  the 
Pacific  Ocean  in  the  vicinity  of  Station  "P".  This  flow 
consists  cf  the  Best  Hind  Drift  and  Subarctic  Current  which 
then  diverge  tc  form  the  southward  flowing  California  Cur¬ 
rent  and  the  northward  flowing  Alaskan  Current  (Tabata, 


B.  HATER  HASSES  AND  FRONTS 

The  Northeast  Pacific  can  be  divided  into  three  major 
regions  based  upon  their  water  characteristics.  Near  the 
vicinity  cf  CKS-P  and  northward  is  the  Subarctic  Pacific 
Hater  Mass  in  which  the  upper  10  0-150  o  layer  is  composed  of 
cold  and  lew-salinity  water.  Tc  the  south  warm  and  highly 
saline  water  exists  in  the  upper  layers  of  the  Subtropic 
Hater  Mass.  In  between  and  along  the  coast  lies  a  region  in 
which  characteristics  frem  both  water  masses  are  combined. 
This  is  known  as  the  Transition  Zone  where  there  is  a  marked 
gradation  of  oceanographic  properties  from  one  region  tc 
another,  with  the  boundaries  separating  these  water  masses 
having  complicated  thermohaline  structures  (Roden,  1970)  . 
The  lccaticn  cf  the  Subarctic  Front  between  the  two  water 
masses  remains  fairly  constant,  from  40N  to  45N,  and  is 
largely  dependent  on  the  prevailing  winds  (Roden,  1975). 

1.  SUBARCTIC  HATER  MASS  CHARACTERISTICS 

The  tain  feature  of  the  vertical  structure  of 
oceanic  properties  is  the  existence  of  an  upper  zone,  a 
halccline  and  a  lower  zone  (Oda,  1963)  .  The  temperature, 
salinity  and  density  profiles  of  these  three  zones  for  sum¬ 
mer  and  winter  are  depicted  in  Fig.  2. 


WINTER  REMNANT 
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The  upper  zone,  which  extends  frcm  the  surface  tc  a 
depth  of  about  ICO  m,  generally  has  a  iow  salinity  level  of 
less  than  33  g/kg  with  the  lowest  values  being  closest  tc 
the  coastal  regions  in  the  late  spring  ((Ida,  1963)  .  In  the 
winter,  the  water  in  this  zone  is  fhorougnly  mixed  with  all 
physical  properties  becoming  homogeneous.  During  the  sum¬ 
mer,  due  tc  lighter  surface  winds,  a  shallow  ischaline  and 
isothermal  layer  exists  from  10  to  30  m  below  the  sea  sur¬ 
face.  Below  this  mixed  layer  is  a  strong  negative  thermo- 
cline  with  tenperature  decreasing  as  much  as  8C  in  20  m;  in 
the  same  zona  a  small  halccline  extends  above  another  neaz- 
isckalire  layer  to  the  bottom  of  the  upper  zone  (Tabata, 
1965)  . 

Below  the  upper  zone  lies  the  main  haloclina  which 
is  permanent,  as  the  salinity  increases  by  as  much  as  1  g/kg 
from  100  tc  200  m  (Tabata,  1960)  .  The  temperature  generally 
decreases  with  depth  in  this  region,  but  at  times  a 
temperature  .nvsrsion  occurs  to  a  depth  of  about  150  m. 
This  positive  temperature  gradient  is  in  a  stable  layer 
(density  increases  with  depth)  since  it  is  located  within 
the  main  halocline  (Uda,  1963). 


In  the  lower  zone,  below  the  nalocline. 


the  proper¬ 


ties  of  the  water  change  gradually  with  depth.  Salinity 
increases  slightly  by  about  0.6  g/kg  in  1000  ta  while  the 
temperature  decreases  to  a  value  of  approximately  2.8C  at  a 
depth  of  1000  a  (Tabata,  1S65). 

The  variations  of  salinity  are  most  pronounced  in 
the  upper  zone  and  halocline  and  are  less  in  the  lower  zone. 
The  annual  salinity  variations  in  the  upper  zone  are  charac¬ 
terized  by  a  decrease  in  the  upper  30  m  in  the  late  summer 
or  early  fall  with  a  maximum  in  late  winter  or  early  spring. 
These  seasonal  fluctuations  result  primarily  from  the  evapo¬ 
ration  and  precipitation  imbalances  common  to  the  region. 
Annual  variation  is  similar  at  great  depths,  but  lags  the 
surface  layers  in  time  (Tatata,  1965)  . 

At  the  surface,  several  features  are  observed  in  the 
salinity  pattern.  Basically,  the  trend  of  the  surface 
isohalines  appears  to  be  east-west  throughout  the  year  with 
major  changes  occurring  along  the  coast.  Also,  a  permanent 
salinity  maximum  of  about  33  g/kg  exists  near  the  center  of 
the  Gulf  cf  Alaska,  south  of  Kodiak  Island.  The  salinity 
decreases  southward  to  a  minimum  of  about  32.6  g/kg  in  the 
vicinity  of  SON  and  then  increases  further  south  into  the 


Subtropical  Mater  Mass  (Uda,  196  3) 


The  "three-zone"  structure  of  the  water  coin  nr  is 
characteristic  of  the  deep  ocean  region  of  the  Subarctic 
Domain.  However,  it  does  not  always  exist  in  near-coastal 
waters  as  current  systems,  coastal  winds  and  fresh  water 
runoff  from  land  can  modify  the  profiles.  As  shown  in  Pigs. 
3  and  4,  this  coastal  regime  primarily  alters  the  upper 
zone.  Along  the  western  United  States  this  domain  exists 
out  tc  about  130H.  The  lower  zone  also  becomes  affected  by 
the  proximity  to  land  of  the  California  Undercurrent 
(Dcdimead  et  al. ,  1962). 

2.  SUETBCEIC  MATES  MASS  CHA  B  ACTER I5TICS 

To  the  south  of  the  Subarctic  region  is  the  Sub¬ 
tropic  Mater  Hass  which  is  basically  warmer  and  more  saline 
than  the  Subarctic  Mater.  In  the  upper  layer  the  high  val¬ 
ues  cf  salinity  (greater  than  34  g/kg)  due  to  higher  evapo¬ 
ration  decrease  to  a  minimum  at  200  to  600  m  and  are  typical 
of  the  values  throughout  the  year.  Sea  surface  temperatures 
remain  fairly  constant  with  an  isothermal  upper  zone  averag¬ 
ing  about  20  m  in  thickness  below  the  sea  surface,  due  tc 
the  prevailing  winds.  Eelow  this  zone  exists  a  negative 
permanent  thermocline  reaching  down  to  a  depth  of  500-600  m. 
Also,  there  is  no  characteristic  halociine  or  major 
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Figure  3.  Upper  Zone  Domains  (From  Dodimeai  e 


Zone  Domains  (From  Dodimead 


temperature  inversion  as  there  was  in  the  water  mass  to  the 
north  (Tatata,  1965).  The  deep  water  structure  is  basically 
continuous  and  similar  from  the  Subarctic  to  Subtropic 
regions  (Oda,  1963)  . 

3.  transition  zone 

The  Transition  Zone,  located  in  the  central  and 
eastern  North  Eacific  between  45N  and  35N,  separates  Subarc¬ 
tic  Hater  frcm  Subtropic  Hater.  It  consists  of  the  mixed 
water  originating  in  the  warm,  saline  Kurcshio  and  the  cold, 
low  salinity  Cyashio  currents  as  they  move  eastward  in  the 
West  Wind  Drift  (Oda,  1963). 

The  width  of  this  zonal  belt,  which  is  about  2  tc  4 
degrees  cf  latitude,  and  the  strength  of  its  boundaries,  are 
primarily  determined  by  the  wind  stress  at  the  sea  surface. 
The  northern  boundary  is  characterized  by  many  temperature 
inversions  ar.d  by  the  gradual  disappearance  of  the  Subarctic 
halocline.  At  the  southern  boundary  the  vertical  structure 
becomes  basically  isohaline  with  a  strong  negative  thermo- 
cline  in  the  upper  layer  (Seder.,  1970)  . 

The  Transition  Zone  is  primarily  a  property  cf  the 
upper  ocean,  as  its  characteristics  fade  below  a  few  hundred 


meters 


At  the  surface  this  region  has  a  distinctive 


meridional  (N-S)  temperature  gradient  of  approximately  1.GC 
per  1001cm  as  it  stretches  eastward  across  the  Pacific  Ocean. 
The  horizontal  surface  salinity  pattern  exhibits  similar 
structure  alcng  a  latitude  band  throughout  the  year  with 
values  increasing  as  much  as  0.8  g/lcg  meridionally  (Uda, 
1963)  . 

The  upper  layers  are  often  characterized  by 
temperature  inversions  which  can  occur  because  of  the  exis¬ 
tence  of  the  associated  halocline.  These  inversions 
decrease  in  magnitude  and  slowly  increase  in  depth  toward 
lower  latitudes. 

The  vertical  distribution  of  properties  at  interme¬ 
diate  depths  shows  that  the  salinity  becomes  effectively 
isohaline,  with  a  slight  minimum  developing  below  300  m  at 
the  southern  boundary  (Uda,  1  963). 

4.  DENSITY  STRUCTURE  AND  THE  SUBARCTIC  FRONT 

The  density  cf  sea  water  is  calculated  primarily 
from  its  temperature  and  salinity  at  depth.  In  the  Subarc¬ 
tic  waters,  salinity  governs  the  density  structure.  The 
seasonal  variability  illustrates  some  of  tae  features  that 
occur  in  the  vertical  profiles.  In  winter,  basically  iso- 
pycnal  conditions  are  found  throughout  the  upper  zone  as 
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temperature  and  salinity  are  well  mixed.  In  summer,  a  minor 
pycnccline  corresponding  to  the  shallow  thermocline  is 
observed.  However,  a  large  and  permanent  pycnoclir.e  occurs 
along  with  the  main  halocline.  In  the  lower  zone,  the  den¬ 
sity  increases  gradually  with  depth  (Dda,  1963) . 

In  Subtropic  Hater,  the  density  is  controlled  by  the 
temperature.  Here  the  warm  surface  layer  of  about  10  to 
30  m  is  isopycnal  from  the  mixing  of  the  prevailing  winds. 
A  permanent  pycnccline  exists  at  lower  depths  due  to  the 
thermocline  in  the  region  (Dda,  1963)  . 

The  Transition  Region  has  a  complex  density  struc¬ 
ture  as  numerous  minor  density  fronts  occur  at  depth  due  to 
the  temperature  and  salinity  inversions. 

The  Subarctic  Front  in  the  eastern  North  Pacific  is 
a  narrow,  meandering  band  located  in  the  latitude  range  of 
40 N  tc  45N.  This  front  has  a  unique  structure  in  its  lack 
of  a  density  front  in  the  upper  100  m.  The  strong  horizon¬ 
tal  temperature  and  salinity  gradients  are  in  almost  com¬ 
plete  balance.  Below  this  upper  layer  a  density  front  of 
moderate  intensity  slopes  slightly  to  the  south  (Roden, 
1975)  . 
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The  mixed  layer  depth  can  strongly  change  across  the 
Subarctic  Frcrt.  Tc  the  north,  the  mixed  layer  extends  to 
the  top  of  the  halocline  at  about  100  m.  To  the  south, 
there  is  nc  halocline  and,  during  the  stormy,  winter  season, 
the  mixed  layer  depth  can  extend  to  almost  300  m.  In  the 
summer  months  the  seasonal  thermocline  controls  the  depth  of 
the  mixed  layer  and  there  is  little  difference  across  the 
Subarctic  Front  (Eoden,  1S75). 


III.  THERMAL  STRUCT USE 


A.  THERMAL  STRUCTURE  IN  THE  NORTHEAST  PACIFIC  OCEAN 

The  interaction  between  the  sea  and  atmosphere  is  the 
dominant  process  determining  temperature  structure  in  the 
upper  layers.  Below,  temperature  is  a  fairly  constant  prop¬ 
erty  with  the  structure  being  determined  by  internal  pro¬ 
cesses,  such  as  advection  and  mixing. 

Two  heating  cycles  occur  in  the  ocean:  diurnal  and  sea¬ 
sonal.  In  the  daily  cycle  the  surface  layer  is  heated  dur¬ 
ing  the  daylight  hours  and  cooled  during  the  night.  In  the 
higher  latitudes  during  the  heating  season,  i.e.,  the  middle 
six  months  of  the  year,  a  daily  net  heat  gain  occurs  (Oda, 
1963).  Fig.  5  illustrates  the  annual  cycle  of  heating  and 
coding  and  its  relation  to  the  thermal  structure.  In 
March,  at  the  end  of  the  coding  season,  the  waters  are 
nearly  isothermal  to  the  top  of  the  haiocline  located  at 
about  100  m  in  depth.  With  the  beginning  of  the  heating 
seascn  in  April,  the  heat  is  accumulated  as  small 
temperature  inversions  are  gradually  mixed  downward  by  light 
winds.  During  subsequent  short  periods  of  high  winds  these 


transients  are  quickly  eroded  and  a  shallow  isothermal  sur¬ 
face  layer  results. 

Eelow  the  surface  layer,  a  sharp  negative  chermocline  is 
formed  between  30  and  60  m  in  a  temperature  transition  layer 
adjacent  to  the  deeper  cold  waters.  Prom  the  and  cf  the 
heating  season  until  the  waters  are  again  nearly  isothermal 
in  March,  the  thermocline  deepens  to  about  100  m. 

In  this  ocean  region,  below  the  thermocline  there  is 
usually  a  coder  layer  nearly  120  m  in  depth  which  is  sta¬ 
ble,  since  it  occurs  within  the  main  halociine.  Below  the 
cooler  layer  there  may  be  an  increase  in  temperature  with 
depth,  an  inversion.  Temperature  inversions  can  be  classi¬ 
fied  into  three  types  according  to  their  method  of  formation 
and  behavior.  The  first  type,  which  is  the  most  common  in 
the  eastern  North  Pacific,  is  a  result  of  the  seasonal  cycle 
of  heating  and  cooling.  The  upper  layers  in  an  extremely 
cold  winter  can  become  colder  than  the  waters  in  the  halo¬ 
ciine  below  and  a  positive  temperature  gradient  is  formed. 
With  the  advent  of  the  following  heating  season  the 
temperature  cf  the  surface  waters  increase  causing  a 
(negative)  thermocline  tc  form  above  the  positive  thermal 
gradient.  The  temperature  minimum  thus  generally  is  located 


just  above  tie  haiocline  at  120  ra  and  changes  very  little 
daring  the  course  of  the  summer  heating  period.  During  the 
following  winter  season  the  upper  thereoclins  decays  and  the 
temperature  inversion  may  be  erased  (Uda,  1963). 

Another  process  which  can  cause  temperature  inversions 
is  the  formation  of  sea-ice  in  the  western  Pacific.  The 
remaining  cold  (-1C)  ,  high  salinity  surface  layer  sinks  tc  a 
level  of  stability  in  the  existing  haiocline  and  is  carried 
eastward  in  the  West  Wind  Drift  where  it  completely  mixes 
with  other  waters  before  reaching  140W.  Depending  cn  the 
amount  of  heating  of  the  surface  layers  during  the  transit, 
the  temperature  inversion  may  or  may  not  persist  throughout 
the  year  (Eoden,  1964). 

A  third  process  results  in  a  subsurface  temperature  max¬ 
imum  when  the  warm,  more  saline  Subtropical  Water  intrudes 
under  the  colder  Subarctic  Waters.  The  advecting  water  then 
is  slightly  ceded  and  sinks;  it  remains  warmer  than  the 
surroundings,  however,  and  forms  the  temperature  maximum  at 
depth.  Its  upper  boundary  is  generally  limited  by  the  sta¬ 
bility  criterion  of  the  main  haiocline  (Uda,  1963)  . 

As  the  ceding  period  begins  in  late  September,  the  sur¬ 
face  waters  are  cooled  and  mixed  downward  by  convection  as 
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well  as  by  the  strong  winds  of  the  fall  and  winter  months. 
The  seasonal  theraocline  sinks  and  narrows  until  it  reaches 
the  halocline.  By  late  winter  the  upper  zone  has  reached 
its  minimum  temperature  of  the  year. 

At  depths  telow  500  a  temperature  variations  are  small, 
seldom  exceeding  standard  deviations  of  0.04C  (Tabata, 
1960).  Fig.  6  shows  the  typical  temperature  structure  at 
Ocean  Station  "P"  for  both  summer  and  winter  months. 

B.  INTERNAL  WAVE  EFFECTS  CN  THERMAL  STRUCTURE 

In  many  vertical  temperature  structure  analyses,  the 
depths  of  the  isotherms  fluctuate  considerably.  Periodicity 
and  amplitude  of  the  fluctuations  appear  non-systemat ic ,  but 
result  from  the  superposition  of  a  large  number  of  infernal 
wave  trains.  The  limiting  frequencies  range  from  the 
Br unt-Vaisala  ic  inertial;  in  the  Subarctic  Region,  periods 
from  5  minutes  to  at  least  half  a  day  have  been  observed. 
The  vertical  displacement  of  isotherms  is  of  the  order  of 
5  m  at  the  top  of  the  thermocline  and  perhaps  as  large  as 
25  m  at  the  level  of  the  principal  halocline  (Tully  and 
Giovando,  1963) . 
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C.  HORIZONTAL  SURFACE  THERMAL  STRUCTURE 

Hqar  the  center  of  the  Gulf  of  Alaska  is  located  an  area 
of  temperature  minimum  that  coincides  with  the  dynamic  cen¬ 
ter  cf  the  Alaskan  Gyre.  In  summer  it  roughly  overlies  the 
canter  of  the  gyre  as  determined  by  the  salinity  maximum, 
but  in  winter  the  two  are  separated  (Uda,  1963). 

South  from  the  cold  core  dome  to  approximately  35N  the 
therial  field  is  oriented  zonally  except  near  the  ccastal 
regions.  This  configuration  continues  year  round  with  a 
seasonal  range  cf  7.0C  throughout  the  Transition  Region 
(Uda,  1963).  The  horizontal  thermal  gradient  m  the  north- 
south  direction  is  strongest  in  the  Transition-  zone  where 
gradients  cf  IC/IOOkm  are  common,  as  shown  in  Fig.  7. 
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IV.  EXPERIMENT 


A.  OBJECTIVE 

For  ocean  monitoring  from  space  to  become  effective,  the 
interior  metiers  of  the  oceans  must  be  inferred  from  the 
resulting  signature  of  surface  properries.  The  ability  tc 
define  the  subsurface  structure  accurately  by  this  method  is 
becoming  increasingly  utilized  and  understood. 

legeckis  and  Gordon  (1962)  described  the  ability  of  sat¬ 
ellite  infrared  imagery  tc  locate  sea  surface  temperature 
fronts,  and  the  relation  of  these  surface  features  to  sub¬ 
surface  temperature  structure.  For  the  Brazil  Current  and 
two  warm  cere  eddies,  they  noted  the  mixed  layer  depth  cor¬ 
related  with  the  surface  temperature  patterns. 

Bernstein,  Ereaker  and  Whritner  (1977)  examined  the  eddy 
formations  in  the  California  Current  from  combined  ship,  air 
and  satellite  platforms  and  found  that  the  sea  surface 
temperature  patterns  conformed  to  subsurface  warm  and  ccld 
mesoscale  variations.  They  concluded  that  in  unstable 
regions  such  as  the  California  Current,  surface  and  main 
thermccline  temperature  distributions  have  relationships  to 
each  ether.  Thus,  in  such  a  region,  remotely  sensed  sea 


surface  temperature  patterrs  would  be  able  to  interpret  sub¬ 
surface  circulation. 

From  these  recent  examples  of  satellite-determined  sub¬ 
surface  structure  and  frcm  other  oceanographic  literature, 
it  appears  that  the  regicr.s  of  the  world's  oceans  where 
large  instabilities  occur  (i. e. ,  boundary  currents,  coastal 
upwelling  zcnes,  etc.),  the  subsurface  dynamic  features  ara 
readily  deduced  from  the  surface  observations.  Areas  of  the 
seas  which  are  quiescent  and  usually  located  in  the  central 
oceans  have  net  been  as  easy  to  describe  accurately.  In 
these  regions,  the  use  cf  either  dynamical  or  empirical 
methods  cculd  be  the  only  means  of  forecasting  the  thermal 
structure  cf  the  cceans. 

Numerical  models  are  based  cn  dynamical  principles  and 
will  ultimately  provide  the  basis  for  predicting  oceanic 
processes.  Empirical  methods  have  been  determined  in  the 
past  for  certain  regions  and  for  limited  time  periods,  tut 
do  net  appear  tc  be  valid  universally.  If  relationships 
could  be  devised  for  certain  water  masses  or  oceanic 
regimes,  then  the  possibility  exists  of  piecing  together  the 
subsurface  thermal  features  in  the  different  regions  by  uti¬ 
lizing  remote  sensing  to  define  certain  surface 
measurements . 


Toward  this  goal,  the  first  objective  of  this  thesis  is 
to  determine  what  particular  relationships  exist  between  sea 
surface  temperature  and  subsurface  thermal  structure  for  a 
relatively  "quiet”  region  in  the  northeast  Pacific  Ocean. 
By  examining  the  climatology,  an  empirical  set  of  equations 
could  then  be  derived  between  sea  surface  temperature,  mixed 
layer  depth,  thermocline  gradient  and  the  depths  of  promi¬ 
nent  points  within  the  structure. 

The  second  objective  is  to  produce  a  predicted  vertical 
temperature  prefile  from  the  best  fit  regression  equations 
for  a  particular  water  region  and  SST  at  a  future  time  and 
different  location.  when  these  profiles  are  compared  to 
actual  temperature-depth  profiles,  accuracy  of  this  statis¬ 
tical  method  will  be  determined. 

B.  E AT A  ACQUISITION 

Curing  the  period  5  tc  14  September  1977,  the  U.S.  Naval 
oceanographic  survey  vessel  USNS  SILAS  SENT  conducted  ocean¬ 
ographic  surveys  to  determine  the  thermohaline  and  sound 
velocity  structures  across  the  Califiornia  Current,  the 
Transitior  Zcre  and  Subtropical  Mater.  For  this  thesis,  the 
only  data  utilized  were  along  a  1155  km  meridional  track 
south  cf  Cesar.  Weather  Station  ”P",  as  displayed  in  Fig.  8. 
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The  data  consisted  of  28  vertical  temperature  profiles 
spaced  approximately  42  km  apart;  they  spannad  the  period  9 
to  11  September  1977. 

Heather  conditions  during  this  period  were  governed  by  a 
weak  stationary  low  pressure  system  with  mostly  low  cloud 
cover  and  periods  of  light  drizzle  and  showers.  The  wind 
speed  varied  from  5  to  20  knots  from  several  directions 
causing  wave  heights  of  approximately  3  to  6  m  along  the 
entire  track  (FNOC,  1977). 

Expendable  bathythermographs  (XBT's)  were  used  to  define 
the  thermal  structure.  The  resulting  temperature  traces 
plotted  onboard  ship  were  coded  or.  bathythermograph  leg 
sheets,  which  were  utilized  as  the  data  points  for  recon¬ 
struction  cf  the  thermal  profiles.  Characteristic  of  this 
procedure  is  the  loss  of  the  fine  thermal  structure,  as  only 
prominent  points  cn  the  temperature  profile  are  recorded  to 
give  the  general  shape  and  location  of  important  features. 

Another  prchlem  associated  with  this  method  of  data 
acquisition  is  the  varying  probability  that  the  observer  on 
the  vessel  at  the  time  of  the  X BT  drop  correctly  chose  the 
proper  points  tc  reflect  adequately  the  mixed  layer  depth, 
thermccline  gradient  and  temperature  inversions.  Incorrect 


recording  cf  the  values  alcng  with  malfunctioning  in  5  no 
10*  cf  the  XBTs  can  also  fce  a  source  of  potential  error. 

C.  ASSESSMENT  CF  TEMPERATURE  STRUCTURE 

Tc  begin  the  analysis  cf  subsurface  thermal  structure, 
realizing  the  "crude’'  nature  of  the  data,  the  temperature- 
depth  points  were  plotted  cn  a  computer  graphics  system  and 
the  28  thermal  profiles  were  reconstructed.  From  these  pro¬ 
files  a  vertical  cross-section  of  the  meridional  temperature 
structure  was  drawn.  Fig.  5.  It  was  noted  immediately  that 
numerous  temperature  inversions  existed  alcng  the  northern 
and  central  regions  of  the  track,  but  they  ceased  at  approx¬ 
imately  4  IN.  This  subsurface  border  clearly  defines  the 
Subarctic  Frcnt  which  divides  the  Subarctic  from  Subtropic 
Waters;  while  at  the  surface  there  is  no  indication  cf  the 
division.  The  temperature  inversions  varied  in  strength, 
but  generally  shewed  an  increase  in  depth  in  a  southward 
direction. 

The  upper  mixed  layer  displayed  little  variablility  with 
few  transients  fo.-nd  within  the  structure.  There  were  no 
large  eddies  rcted  in  the  surface  layer  and  the  mixed  layer 
depth  appeared  relatively  constant  along  the  entire  track. 
Weak  mescscale  features  are  probably  present,  but  the 
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spatial  resolution  prohibits  their  detection 


The  t  her  de¬ 


cline  was  also  fairly  constant  in  depth  with  the  strongest 
gradients  between  30  to  6C  m.  The  isotherms  in  the  lower 
layers  characteristically  deepened  to  the  south,  with  sev¬ 
eral  areas  of  large  variability  located  near  the  temperature 
inversions.  The  large  warm  intrusion  located  at  47N  may  be 
caused  by  errors  in  the  data,  as  a  structure  this  size  in 
the  Subarctic  Region  does  not  appear  to  be  congruent  with 
the  ether  features,  and  is  not  mentioned  in  the  research 
li terature . 

To  assess  further  the  vertical  thermal  structure,  promi¬ 
nent  points,  gradients  and  zones  were  defined  on  a  typical 
temperature  profile  from  the  test  region.  Fig.  10  illus¬ 
trates  the  variables  that  were  used  to  analyze  statistically 
the  thermal  profile.  Descriptions  of  these  variables  are: 


•  SST  -  Sea  surface  temperature  was  defined  as  the  sur¬ 
face  temperature  reported  on  the  bathythermograph  logs. 
Since  this  is  the  temperature  recorded  on  the  XbT  trace 
when  the  XET  Drobe  entered  the  water,  a  slight  dispar¬ 
ity  from  remotely-sensed  SST  may  occur  due  to  the 
response  time  of  the  XBT  system. 


•  OLD  -  nixed  layer  depth  was  visually  determined  as  the 
greatest  depth  where  the  isothermal  or  slightly 
negative  temperature  gradient  in  the  surface  layer 
changed  to  the  stronger  gradient  of  the  underlying 
theracclme. 


•  BOT  -  Ecttcm  of  seasonal  thermociine  was  visually 
determined  as  the  depth  where  the  seasonal  thermocline 
markedly  changes  into  a  weaker  negative  temperature 
gradient  fceiew. 
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•  TBO I  -  Temperature  at  the  bottom  of  the  seasonal  zher 
mocline. 


DT  -  Difference  in 
at  the  mixed  layer 
ther  mcdine . 


temperature  between  the  temperature 
deptn  and  at  the  bottom  of  seasonal 


DZ  -  Difference  in  depth  between  the  mixed  layer  depth 
and  bottom  of  seasonal  thermocline. 


•  EOZ  -  Depth  of  the  bottom  of  the  main  thermocline  was 
visually  determined  as  the  depth  where  the  thermocline 
changes  to  an  isothermal  or  positive  temperature 
gradient  *ust  above  the  temperature  inversion  region. 


«  TBUZ  -  Temperature  at  the  bottom  of  the  main  thermo¬ 
cline  . 


•  DTBUZ  -  Difference  in  temperature  between  the 
temperature  at  the  mixed  layer  depth  and  at  the  bottom 
of  main  thermocline. 


•  DZB0Z  -  Difference  in  depth  between  the  mixed  layer 
depth  and  the  bottom  cf  main  thermocline. 


•  T460  -  Temperature  at  reference  deptn  of  460  m  which 

was  selected  arbitrarily  at  maximum  depth  of  XBT 
reports. 


•  DT5  -  Temperature  change  in  upper  5  m  of  thermocline 
below  the  mixed  layer  depth. 


«  DT15  -  Temperature  channge  in  upper  15  m  of  thermocline 
below  the  mixed  layer  depth. 


•  DT50  -  Temperature  change  in  upper  50  m  of  thermocline 
below  the  mixed  layer  depth. 


•  DT100  -  Temperature  change  in  upper  100  m  cf  thermc- 

cline  below  the  mixed  layer  depth. 


•  DIS  -  Distance  from  northern  most  XBT  drop. 
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This  assessment  defines  the  general  structure  primarily 
in  the  upper  layers  above  the  permanent  halocline.  The 

detail  of  the  smaller  structure  such  as  transients  in  the 
mixed  layer  and  multiple  ther moclinss  are  not  included 
because  they  are  variable  and  minor  with  respect  to  the 
principal  structure.  The  temperature  inversion  feature 
below  the  main  thermocline  was  not  characterized  as  its 
occurrence  in  the  region  varies  in  location  and  intensity 
during  the  test  period.  In  those  zones  within  the  Subarctic 
Region  where  a  sub-t hermccline  temperature  maximum  exists 
permanently  during  the  heating  season,  additional  variables 
should  be  used  tc  define  this  part  of  the  thermal  profile. 

Initially,  the  sea  surface  temperature  gradient  was  used 
as  a  variable  in  the  analysis.  However,  problems  occurred 
ir.  computing  consistent  gradients  along  the  track  as  the 
calculated  values  were  unrealistic  and  inaccurate  compared 
with  actual  conditions  at  the  time;  this  variable  was 
dropped  from  this  study. 

These  verities  provided  the  statistical  base  from  which 
all  correlations  and  regression  equations  were  derived.  The 
values  could  then  distinguish  the  fluctuating  and  transient 
features  from  the  median  thermal  structure.  Tables  I,  II 
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and  III  list  the  values  from  the 


by  XET  drop  number. 


TABLE  I 

locations  and  Times  of 


XBT 

C  ATE^TIME  i 

(Z) 

LAT 

<N) 

339 

0906  13 

SEP 

77 

47. 

98 

340 

0SC804 

SEP 

77 

47. 

65 

341 

C91003 

SEP 

77 

47. 

20 

342 

C  5  1200 

SEP 

77 

46. 

82 

343 

091800 

SEF 

77 

46. 

07 

244 

0S2000 

SEP 

77 

45. 

70 

345 

092205 

SEP 

77 

45. 

33 

346 

100000 

SEP 

77 

45. 

02 

347 

1C0308 

SEP 

77 

44. 

73 

248 

1005C8 

SEP 

77 

44. 

37 

349 

100708 

SEP 

77 

44. 

00 

3  50 

1CC904 

SEP 

77 

43. 

63 

351 

1 0 1 1  CO 

SEP 

77 

43. 

28 

352 

1 C 1 400 

SEP 

77 

43. 

07 

353 

101600 

SEP 

77 

42. 

70 

354 

1C  1800 

SEP 

77 

42. 

32 

355 

102000 

SEP 

77 

41. 

93 

356 

1 1C002 

SEP 

77 

41. 

55 

357 

110203 

SEP 

77 

41. 

17 

358 

1104C6 

SEP 

77 

40. 

75 

359 

1 10604 

SEF 

77 

40. 

37 

360 

111000 

SEP 

77 

39. 

97 

361 

1 1 1200 

SEP 

77 

39. 

60 

362 

111400 

SEP 

77 

39. 

22 

363 

111600 

SEP 

77 

38. 

83 

364 

1  12304 

SEP 

77 

38. 

30 

365 

120107 

SEP 

77 

37. 

93 

366 

1203C7 

SEP 

7  7 

37. 

58 

analyzed  thermal  profiles 


Measurements 

LONG  (») 

DISTANC 

144. 98 

0 

144. 98 

36 

144.98 

86 

144. 82 

128 

145. 05 

212 

145. 04 

253 

145.04 

294 

145. 04 

329 

145. 00 

361 

144. 98 

401 

144. 98 

442 

144. 98 

483 

144.97 

522 

144. 95 

546 

144. 93 

587 

144. 95 

629 

144. 98 

672 

145.00 

714 

145. 00 

757 

145.00 

803 

145. 00 

846 

144. 98 

890 

145.  10 

931 

145. 25 

973 

145. 37 

1017 

145. 47 

1076 

145. 53 

1117 

145. 67 

1156 

TABLE  II 

Values  of  SSI,  MID  r  Thermoclir.a  Gradianrs  and  T460 


X3T 

SST  (C) 

KID  (m) 

DTS (c) 

£115  (c) 

DT 50  (c ) 

DTI 00  {£) 

T4 6 0  (c) 

3  39 

13.4 

_  “  -i 

1.5 

4.6 

7.1 

6.6 

4.0 

340 

13.6 

-24 

0.3 

1.  5 

6.6 

7.4 

4.1 

341 

14.  1 

-25 

2.5 

4.3 

6.0 

6.  2 

5.0 

342 

14.4 

-25 

2.3 

5.6 

7.1 

7.6 

4.2 

343 

15.  1 

-19 

2.4 

4.8 

7.3 

8.  1 

4.6 

344 

15.3 

-32 

0.  3 

4.6 

7.7 

7.4 

4.9 

3  45 

15.6 

-29 

3.  1 

4.9 

7.3 

7.  9 

6.2 

3  46 

15.7 

-19 

0.  4 

4.7 

6.7 

7.9 

4.7 

3  47 

16.6 

-27 

0.2 

4.  8 

7.9 

9.  1 

4.7 

348 

16.2 

-20 

0.  5 

4.  7 

7.4 

8.  4 

4.8 

349 

17.0 

-21 

0.4 

5.3 

8.0 

8.  5 

4.9 

350 

18.  1 

-20 

0.  2 

6.2 

9.1 

10.  3 

4.9 

351 

18.3 

-30 

3.  6 

4.2 

7.8 

10.  2 

5.0 

352 

18.  1 

-20 

0.3 

5.8 

8.4 

8.7 

5.1 

353 

18.4 

-20 

0.7 

4.0 

9.0 

9.9 

5.0 

354 

18.1 

-20 

0.6 

5.4 

8.6 

9.3 

5.2 

355 

19.3 

-21 

0.  6 

5.4 

9.3 

10.  3 

5.3 

3  56 

19.7 

-15 

0.2 

6.9 

9.6 

10.8 

5.3 

357 

19.5 

-23 

0.6 

6.4 

9.7 

10.  7 

5.3 

358 

20. 1 

-16 

0.3 

2.  1 

9.4 

10.3 

5.6 

359 

20.3 

-22 

0.8 

7.0 

10.0 

10.8 

5.6 

360 

20. 1 

-16 

0.5 

2.6 

8.3 

10.3 

5.9 

361 

20.5 

-30 

2.4 

6.1 

9.3 

10.4 

5.5 

362 

21.5 

-22 

0.3 

5.7 

10.5 

11.6 

5.8 

363 

21.4 

-32 

3.8 

6.4 

9.8 

11.4 

5.8 

364 

22.5 

-27 

0.4 

6.2 

9.8 

11.1 

6.5 

365 

22.3 

-22 

0.4 

5.6 

7.6 

9.6 

6.9 

366 

22.7 

-20 

3.  3 

5.2 

8.  1 

10.  4 

6.8 
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TABLE 

III 

Values  of  Th 

srmocline  Temperat 

urss  and 

Dept  hs 

XBT 

EOT (m) 

TECT  (C) 

DT  (c) 

CZ  (m) 

BOZ  (m) 

TBOZ  (c) 

DTBUZ(c) 

DZEUZ  (m ) 

339 

-88 

5.4 

7.4 

55 

-95 

5.3 

7.  2 

62 

340 

-30 

6.5 

6.8 

56 

-109 

5.9 

7.4 

85 

34  1 

-70 

7.8 

5.9 

45 

-110 

7.4 

6.  3 

85 

34  2 

-79 

6.6 

7.2 

54 

-128 

6.2 

7.6 

103 

34  3 

-81 

6.6 

8.0 

62 

-138 

6.5 

8.  1 

1  19 

344 

-81 

7.  1 

7.7 

49 

-112 

6.7 

8.  1 

80 

345 

-68 

7.8 

7.1 

39 

-112 

7.0 

7.  9 

83 

346 

-57 

8.4 

6.8 

38 

-100 

7.3 

7.9 

81 

34  7 

-68 

8.6 

7.8 

41 

-114 

7.  3 

9.  1 

87 

34  8 

-65 

6.4 

7.3 

45 

-143 

7.6 

8.  1 

128 

34  9 

-50 

9.1 

7.5 

29 

-109 

8.  1 

8.  5 

88 

350 

-56 

8.8 

9.0 

36 

-123 

7.  5 

10.3 

103 

35  1 

-57 

10.7 

7.0 

27 

-148 

7.  2 

10.  5 

1  18 

35  2 

-73 

9.0 

8.5 

53 

-158 

7.3 

9.7 

138 

353 

-67 

9.4 

8.8 

47 

-132 

7.7 

10.  5 

112 

354 

-63 

9.3 

8.6 

43 

-120 

8.6 

9.  3 

100 

35  5 

-59 

9.9 

9.1 

38 

-130 

8.6 

10.  4 

109 

35  6 

-54 

9.9 

9.5 

39 

-134 

8.5 

10.  9 

119 

35  7 

-6  2 

10.1 

9.1 

39 

-127 

8.4 

10.  8 

104 

358 

-50 

11.3 

8.6 

34 

-120 

9.  3 

10.6 

1  04 

35  9 

-53 

10.8 

9.4 

31 

-95 

9.4 

10.  8 

73 

360 

-78 

10.5 

9.5 

62 

-180 

9.3 

10.7 

164 

36  1 

-99 

9.5 

10.2 

69 

-127 

9.3 

10.  4 

97 

362 

-68 

11.0 

10.4 

46 

-157 

9.4 

12.0 

135 

363 

-100 

10.2 

1  1  .0 

68 

-155 

9.7 

11.5 

123 

364 

-91 

12.2 

10.0 

64 

-173 

10.5 

11.7 

144 

365 

-100 

12.5 

9.3 

78 

-196 

10. 7 

11.  3 

174 

366 

-100 

12.2 

10.2 

80 

-188 

10.7 

11.7 

168 

P 


r. 

t 


r 
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To  examine  the  variations  and  trends  of  the  individual 
therial  profile  features  along  the  track,  graphs  of  the 
variables  as  related  to  distance  and  latitude  are  plotted  in 
Appendix  A.  These  diagrams  illustrate  some  of  the  concepts 
that  were  discussed  in  previous  chapters  pertaining  to  the 
variability  in  thermal  structure  within  a  particular  region. 

The  sea  surface  temperature  over  the  entire  track  as 
shewn  in  Fig.  A.1  increases  linearly  to  the  south  where  the 
mean  gradient  was  0.8C/100km.  The  strongest  gradient 
occurred  in  the  vicinity  of  44U  and  reached  a  value  of 
2. 3C/100ka.  Fig.  It  shews  the  sea  surface  temperature  in 
the  north  Pacific  as  compiled  from  ship  reports  and  clima¬ 
tology  by  Fleet  numerical  Oceanography  Center  for  10  Septem¬ 
ber  1977.  In  the  region  of  data  acquisition,  the  meridional 
temperature  gradient  and  surface  isotherm  structure  closely 
compare  tc  the  surface  temperature  versus  distance  plot 
(Fig.  A.1).  Satellite  infrared  imagery  during  this  time 
period  was  obstructed  due  tc  the  thick  cloud  cover  over  the 
area;  therefore,  this  irethod  of  sea  surface  temperature 
determination  was  not  used  as  was  desired  (Tabata  and 
Kimber,  1 S79)  . 
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Figure  11 


Pacific  Ocean  Sea  Surface  Temperature  on  10 
September  1977  (Fleet  Numerical  Oceanography 
Center,  1  977). 


Fig.  A.  2  shews  the  depth  of  the  mixed  layer  along  The 
track.  The  mean  depth  was  -22.9  m  with  a  maximum  and  mini¬ 
mum  depth  cf  -32  and  -15  m,  respectively.  The  layer's  gen¬ 
eral  trend  of  a  slight  shallowing  by  about  5  m  in  the 
southward  direction  is  observed  until  40N  is  reached  where 
the  trend  reverses  and  becomes  deeper.  This  location  of 
change  in  the  mixed  layer  depth  trend  accurately  coincides 
with  the  Subarctic  Front. 

The  feur  thermocline  gradient  variables  in  Figs.  A. 3 
through  A. 6  displayed  a  measurable  trend  along  the  track. 
The  variability  of  the  temperature  change  in  the  upper  5  m 
of  the  thermocline  (DT5 )  was  large  in  the  northern  and 
southern  ends  of  the  track.  In  the  center,  with  the  excep¬ 
tion  cf  ens  data  point,  a  mean  value  of  0.3C/5m  was 
observed.  The  boundaries  of  these  three  zones  of  the  DT5 
variable  agreed  with  the  division  of  the  track  into  the 
three  regions,  i.e. ,  Subarctic  Water,  Transitional  Zone  and 
Subtropic  Water.  This  effect  was  noted  in  many  of  the  ther¬ 
mal  structure  variable  plots.  DT15,  DT50  and  DT100  all 
increased  with  decreasing  latitude,  but  the  general  trends 
differed  between  each  variable.  DT15  showed  the  largest 
range  of  variability  with  only  a  slight  increase  in 
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temperature  gradient,  while  DT50  and  DT100  both  increased  at 
a  greater  rate.  Hhen  plotted  collectively,  as  shown  in  Fig. 
12,  the  relationships  become  apparent.  Generally,  with 
increasing  sea  surface  temperature,  the  negative  temperature 
gradient  of  the  thsrmccline  increases  slowly  at  a  relatively 
constant  rate. 

Two  ether  depth-dependent  variables  on  the  vertical 
temperature  profile,  besides  the  mixed  layer  depth,  are  the 
bottem  of  the  seasonal  and  main  thermoclines,  BOT  and  BUZ, 
respectively.  Figs.  A. 7  and  A. 11  display  these  variables  as 
they  vary  with  latitude.  The  variability  of  the  plots  may 
be  due  in  part  tc  internal  waves  with  typical  amplitudes  cf 
15  m  for  the  ECT  and  30  m  for  the  BUZ. 

The  BCT's  general  trend  southward  along  the  track  to 
about  UON  shewed  a  slight  decrease  in  depth  from  -80  to 
-60  m.  Again,  this  location  coincides  with  the  Subarctic 
Front  as  did  the  HLD.  The  BOT  then  increases  sharply  as  the 
trend  reverses  tc  a  depth  of  -90  m. 

For  the  BUZ,  the  trend  is  just  the  opposite  of  the  EOT 
in  the  Subarctic  Hater  and  Transitional  Zone.  The  depth  cf 
the  bottom  of  the  main  thermocline  increases  slowly  until 
reaching  the  Subarctic  Front,  where  the  maximum  variation  in 
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depth  is  approximately  85  m.  In  Subtropic  Water,  the  trend 
of  the  BOZ  increases  in  depth  at  a  greater  rate  along  the 
southern  end  cf  the  track. 

When  the  KID,  BOT  and  EOZ  are  plotted  together,  as  shown 
in  Fig.  13,  additional  relationships  appear.  In  the  Subarc¬ 
tic  and  Transition  Regions,  north  of  40N,  the  thickness  cf 
the  upper  layers  above  the  bottom  of  the  seasonal  thermo- 
cline  (BOT)  decreased  in  a  southward  direction.  With  the 
SLD  trend  at  a  relatively  constant  depth,  the  only  thermal 
layer  to  change  in  magnitude  appreciably  is  the  thermocline. 
This  suggests  that  the  seasonal  thermocline,  with  its 
negative  temperature  gradient,  becomes  stronger  as  the  SST 
increases  and  as  the  thermocline  thickness  decreases. 
Indeed,  the  temperature  difference  (DT)  and  thickness  (DZ) 
of  the  seasonal  thermocline  displayed  in  Figs.  A. 9  and 
A.  10,  respectively,  agree  with  this  suggestion. 

As  previously  mentioned,  the  depth  of  the  bottom  cf  the 
main  thermccline  gradually  deepens  toward  the  south  along 
the  track.  Figs.  A. 13  and  A. 14  show  the  increasing  differ¬ 
ence  in  temperature  and  depth  from  the  mixed  layer  depth, 
DT  BOZ  and  DZBUZ,  respectively.  These  tendencies,  in  con¬ 
junction  with  the  trend  of  the  BOT,  indicate  that  the 
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negative  temperature  gradient  at  the  lower  part  of  the  main 


thermocline  becomes  weaker  and  thicker  toward  the  south  in 
the  Subarctic  and  Transiticnal  Regions. 

For  Subtropical  Mater,  south  of  40N,  the  BOT  has 

reversed  its  trend  at  the  Subarctic  Front  and  increases  in 

depth.  The  BOZ  continues  to  deepen,  but  at  a  greater  rate 

than  to  the  north.  Temperatures  at  these  two  locations  on 
the  thermal  profile,  IBOT  and  TBOZ,  slowly  increase  linearly 
with  decreasing  latitude  at  nearly  the  same  rate  of 
0.5C/100km.  These  values  suggest  that  the  negative 
temperature  gradient  at  the  lower  part  of  the  main  thermo¬ 
cline  continues  to  weaken  and  thicken,  but  at  a  deeper  depth 
than  in  the  northern  waters.  At  the  upper  layers,  the 

parameters  of  the  seasonal  thermocline,  DT  and  DZ,  increase 
with  increasing  SST. 


D.  STATISTIC  II  EROCEEURE 

^  To  predict  the  subsurface  thermal  structure  from  a  given 

a 

sea  surface  temperature,  a  statistical  approach  was  utilized 
to  determine  the  relationships  between  the  thermal  structure 
►I  variables.  The  data  analysis  was  done  using  the  BIOMED  sta¬ 

tistical  computer  programs  on  the  Naval  Postgraduate 


School’s  IBM  3033  computer  system 


Details  of  the  BIOMED 


computer  programs  used  in  the  analysis  are  provided  ir. 
Appendix  E.  Basic  statistics  (mean  and  standard  devia¬ 
tions),  scatter  diagrams,  correlation  and  linear  regression 
analysis  were  performed  on  all  variables. 

As  part  cf  the  regression  analysis,  the  correlation 
coefficients  fcr  all  combinations  of  pairs  of  the  sixteen 
variables  were  computed.  These  simple  correlations  are  an 
estimate  cf  the  strength  cf  the  linear  association  between 
two  variables  and  do  not  imply  the  cause  for  the  relation¬ 
ship  or  which  variables  affect  the  variation  of  the  other. 

Scatter  diagrams  were  used  in  conjunction  with  the 
regression  analysis  to  assist  in  refining  the  correlations 
between  variables.  By  visually  removing  poor  data  points 
that  detracted  from  trends  within  specific  regions,  the 
analysis  cculd  be  improved.  After  determining  the  best  fit 
correlations,  least  square  linear  regression  equations  are 
computed  between  variables  and  also  with  sea  surface 
temperature  as  the  independent  variable.  This  enabled  a 
vertical  temperature  profile  to  be  determined  by  calculating 
the  values  fcr  the  variables  with  the  only  input  being  the 
sea  surface  temperature. 


The  values  cf  the  16  thermal  structure  variables  from 
the  28  temperature  profiles  along  the  track  were  correlated 
for  all  combinations  with  the  results  shown  in  Table  IV. 
Several  strong  relationships  appear  with  the  correlation 
coefficients  cf  a  few  variables  being  greater  than  0.750. 

The  sea  surface  temperature  had  the  greatest  number  of 
strong  correlations  which  greatly  enhances  the  idea  of  pre¬ 
diction  cf  the  thermal  profile  from  the  SST.  The  highest 

correlation  coefficient  (0.993)  occurred  when  the  SST  was 

/ 

compared  with  the  distance  along  the  track.  This  was 
expected  due  to  the  relatively  constant  north-south  surface 
temperature  gradient  in  this  region  of  the  Northeast 
Pacific. 

ether  strong  correlations  with  the  sea  surface 
temperature  were  the  temperatures  at  the  bottom  of  the  sea¬ 
sonal  and  main  thermocline,  TBOT  (0.944)  and  TBUZ  (0.949), 
respectively.  When  plotted  together,  as  shown  in  Fig.  14, 
the  trends  clearly  appear.  The  variables'  similar  relation¬ 
ships  imply  that  the  entire  vertical  temperature  profile 
"leans"  toward  warmer  temperatures  as  the  small  negative 
thermal  gradient  in  the  deeper  layers  increase  with  increas¬ 
ing  sea  surface  temperature.  The  change  in  temperature  in 


TABLE  IV 

Correlation  Coefficients  of  Thermal  Variables 


hLD  MLO 


SST 

CM 

CM 

• 

SST 

DT5 

-.484 

-.023 

DT5 

DTI  5 

-.121 

.335 

.080 

DTI  5 

DT50 

.157 

.770 

-.211 

.485 

or  50 

DTI  00 

.235 

.905 

-.065 

•  393 

.902 

DT100 

DT 

.128 

.335 

-.026 

.466 

.343 

.867 

DT 

DZ 

-.156 

.305 

.280 

.018 

-.082 

.050 

.425 

SOT 

.459 

-.203 

-.406 

-.055 

.124 

.030  - 

.341 

TBOT 

.296 

.944 

-.065 

.251 

.632 

.322 

.638 

3UZ 

-.202  • 

-.632 

-.124 

-.098 

-.239 

-.4q6  - 

.566 

TflUZ 

.267 

.9^9 

-.035' 

•  334 

.630 

.770 

.792 

I/TBUZ 

.217 

.960 

-.068 

.368 

.842 

.961 

.391 

DZBUZ 

.374 

.673 

-.031 

.068 

.246 

.505 

.551 

DIS 

.217 

.993 

-.012 

.365 

.738 

.873 

.685 

T460 

.114 

.531 

.139 

.265 

.435 

.  5f,9 

.644 

MEAN 

-23.2 

18.1 

1.2 

5.0 

3.3 

9.3 

5.5 

STD  DEV  5. 1 

2.8 

1 .2 

1.3 

1.2 

1.5 

1.3 

DZ 

■.949 

BCT 

.150 

-.041 

TBOT 

■  •643 

.5*4 

-.653 

BUZ 

•  355 

-.234 

.932 

-.654 

TEUZ 

.204 

-.114 

.335 

-.636 

.830 

DT3UZ 

.578 

.401 

.665 

-.933 

.661 

.6j4 

DZBUZ 

.356 

-.251 

.931 

-.633 

.'63 

.934 

.634 

DIS 

.405 

-.3 :: 

.835 

-.t66 

.832 

.706 

•  663 

.346 

46,6 

-72.0 

9-3 

-133-3 

£  .1 

9.6 

110.2 

580.9 

14.3 

15.  < 

1.6 

27.4 

1.4 

1 .6 

28.0 

343.6 

the  upper  50  and  100  m  cf  the  thermocline  also  displayed 
strong  direct  relationships  with  the  SST. 

Two  pairs  of  variables  which  exhibited  strong  negative 
correlations,  and  for  similar  reasons,  were  BOT-DZ  (-0.949) 
and  EUZ-DZEUZ  (-0.983).  These  combinations  indicate  that  a 
relationship  exists  between  both  the  30T  and  BUZ  and  the 
MLD.  This  is  due  to  the  relatively  small  variations  in  the 
mixed  layer  depth  when  compared  to  the  large  fluctuations  of 
the  EOT  and  BCZ.  Therefore,  as  the  bottom  of  the  seasonal 
and  main  thermoclines  varied,  the  depth  difference  between 
these  variables  and  the  MID  similarly  changed.  This  could 
provide  an  important  link  in  determining  the  mixed  layer 
depth  if  a  viable  connection  cculd  be  established  between 
the  EOT  or  BUZ  and  SST.  However,  according  to  the  correla¬ 
tion  coefficients  comparing  the  SST  to  MLD  (0.  224),  to  BOT 
(-0.  203)  and  to  BUZ  (-0.682)  ,  the  strength  of  these  rela¬ 
tionships  was  weak  as  shewn  in  Figs.  15,  16,  and  17.  Corre¬ 
lations  of  the  ether  variables  with  the  mixed  layer  depth 
revealed  little  as  to  possible  methods  to  determine  this 
important  thermal  structure  variable. 

Due  tc  the  noticeable  differences  in  thermal  structure 
characteristics  among  the  various  regions  encountered  alcng 


Figure  14 


Ccaparison  of  S ST,  TBOT  and  TBUZ  Versus 
Latitude . 


LEGEND 


-JO 


-so 


-00 


-40-1 


-70-1 


-110 


(tl**,,*/*»,*l*i,l**f*',% 


/ 


A 


»«**«*' 


A 

/  v 


•V 
/  V 
: 


» 


vV 


/ 


% 

:  / 

V# 


•/ 


LEGEND 

EOT 

mimiMiofo 


*• 

« 

1 

*. 


i. 


A 


\s  \  A 
T  r 


40 


44  43  49 

NORTH  LATITUDE  (DEG) 


41 


Figure  16.  Comparison  of  EOT  and  SST  Versus  Latitude. 
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Figure  17.  Comparison  of  EOZ  and  SST  Varsus  Laninud 


tha  track,  the  data  were  next  divided  into  Subarctic,  Tran¬ 
sitional,  Frontal  and  Subtropical  groupings.  These  divi¬ 
sions  were  determined  from  the  previously  discussed  regional 
characteristics  and  thermal  analysis  of  Chapters  II  and  III, 
and  in  retrospect,  they  proved  to  be  an  accurate  reflection 
of  the  boundary  locations  at  the  time  of  data  acquisition. 

The  Subarctic  Water  characteristics  were  traced  from  the 
northern-most  vertical  temperature  profile  on  the  track  to 
XBT  station  number  345  at  45.33N.  The  Transitional  Zone  was 
observed  from  XET  station  346  to  station  357  at  41.17N  fol¬ 
lowed  by  the  Subarctic  Frcnt  between  stations  358  and  359. 
Beginning  at  XBT  station  360  at  (39.97N),  and  continuing 
south  to  the  end  of  the  track,  was  an  ocean  region  contain¬ 
ing  properties  of  the  Subtropic  Water.  Fig.  18  illustrates 
the  divisions  of  the  track  data  points  into  the  various 
regions. 

To  observe  the  separation  of  the  data  by  groupings  bet¬ 
ter,  and  to  refine  farther  the  correlation  coefficients, 
scatter  diagrams  of  several  of  the  thermal  variables  were 
plotted  using  the  BIO  NED  computer  programs.  Since  the  only 
variables  that  correlated  strongly  with  SST  were  other 
temperature-dependent  variables  such  as  TBOT ,  DT,  TBUZ  and 


DTBUZ,  a  closer  inspection  of  the  depth- dependent  variables 
was  desired.  To  ascertain  whether  the  sea  surface 
temperature  can  define  the  depth  of  certain  variables  is 
necessary  in  the  development  of  the  thermal  profile. 
Therefore,  regional  scatter  diagrams  for  the  three  primary 
depth  dependent  variables  (ML D,  BOT,  BUZ)  and  SST  as  com¬ 
pared  to  the  ether  variables  were  plotted  and  analyzed. 

The  scatter  diagrams  for  sea  surface  temperature  when 
compared  to  the  ether  variables  are  listed  in  Appendix  C. 
Due  to  the  high  correlation  coefficients  of  these  thermal 
variables  the  trends  of  the  data  points  were  clearly  visi¬ 
ble.  Obvious  linear  relations  occur  between  SST  and  DT50, 
DT100,  DT,  TECT,  TBUZ  and  DTBOZ  through  the  various  regions. 
DT5  and  DT15  also  displayed  an  element  of  linearity  when 
several  of  the  mere  distant  data  points  were  removed.  The 
computer-derived  correlation  coefficients  for  these  two 
variables  were  excessively  low  due  to  the  remote  points. 
For  the  variable  DT5  shown  in  Fig.  C.2,  the  data  behaved  in 
a  strong  linear  fashion  in  the  Transitional  Region  as  a 
temperature  change  of  0.3C  in  the  upper  5m  of  the  thermc- 
cline  was  the  mean.  All  cf  the  thermociine  gradient  vari¬ 
ables,  DT5,  DT15,  DT50  and  DT100,  could  then  be  modeled  by  a 


linear  equation  related  tc  the  SST,  at  least  in  the  Tran¬ 
sition  Region. 

Two  scatter  diagrams  which  compare  SST  to  BUZ  and  DZ3UZ, 
Figs.  C.  1C  and  C.13,  illustrate  a  perceivable  linear  trend 
between  the  variables  while  the  correlation  coefficients 
reflect  a  moderate  level  of  association,  (-0.687)  and 
(0.678),  respectively.  Again,  with  several  of  the  remote 
data  points  removed,  an  improved  correlation  coefficient 
could  be  computed  which  better  portrays  the  relationship 
with  SST. 

Other  depth-dependent  variables,  for  which  a  visual 
relationship  was  observed  with  SST  in  the  scatter  diagrams, 
were  for  ECT  and  DZ,  Figs.  C.6  and  C.9.  Both  exhibited  a 
"V"  or  wedge  shape  with  data  for  the  Subarctic  and  Subtropi¬ 
cal  Haters  on  either  wing  and  data  for  the  Transitional  Zens 
in  the  center.  If  the  northern  or  southern  regional  data 
were  removed,  then  a  higher  correlation  could  be  obtained 
utilizing  that  from  the  Transitional  Region  alone.  However, 
with  the  visual  trend  in  this  region  closely  following  the 
trend  of  data  from  the  Subarctic  Hater,  data  from  the  north¬ 
ern  end  of  the  track  could  be  used  to  assirt  in  the  accurate 
determination  of  relationships  with  the  SST.  The 
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combination  of  data  from  these  two  regions  is  reasonable 
since  the  location  of  the  Subarctic  Front  is  on  the  southern 
boundary  cf  the  Transition  Zone,  and  the  water  characteris¬ 
tics  are  similar.  Therefore,  the  Subtropic  Water  data 
points  should  be  removed  to  derive  the  correlations  with 
data  from  only  the  two  northern  regions. 

The  two  data  points  in  the  regional  scatter  diagrams  for 
the  Subarctic  Front  consistently  followed  the  Transitional 
Zone  data  trends.  As  shewn  in  Appendix  C,  for  the  thermal 
variables  BOT,  DZ,  BUZ  and  DZBUZ,  the  Front  provided  the 
sharp  break  between  the  Transitional  and  Subtropical 
Regions,  but  in  the  scatter  plots  for  the  remaining  thermal 
variables  the  Subarctic  Front  fell  within  the  trend  of  the 
two  northern  regions. 

The  scatter  plot  comparing  mixed  layer  depth  to  SSI, 
shewn  in  Fig.  C.1,  initially  appears  not  to  reveal  any  pat¬ 
tern  which  could  be  used  to  infer  a  relationship  with  the 
SST.  However,  if  the  Subtropical  Water  data  points  are  not 
included,  a  trend  of  decreasing  depth  is  noticed  and  the 
correlation  coefficient  value  of  0.  224  between  these  two 
variables  should  increase  appreciably. 


The  regional  scatter  diagrams  comparing  mixed  layer 
depth  to  the  ether  variables  did  not  exhibit  any  strong  cor¬ 
relations  and  were  not  included.  However,  several  weak  pat¬ 
terns  were  observed  which  indicated  a  trend  of  the  MLD 
slowly  decreasing  in  depth  in  a  southward  direction  for  the 
northern  two  regions.  All  scatter  plots  displayed  large 
variability  among  tne  data  points  which  agrees  with  the  low 
correlation  coefficients  for  the  variable. 

Appendices  D  and  E  contain  the  regional  scatter  diagrams 
for  ECT  and  EOZ,  respectively.  These  plots  show  a  greater 
linear  dependence  among  the  variables,  and  therefore  higher 
correlation  coefficients,  than  was  shown  by  the  mixed  layer 
depth.  It  most  cases  when  the  Subtropical  Water  data  are 
removed,  leaving  only  data  from  the  Subarctic,  Transition 
and  Frontal  regions,  a  stronger  relationship  was  observed. 

Two  noteworthy  scatter  diagrams  which  had  strong 
negative  correlations  were  B.OT-DZ  and  BUZ-DZBUZ  as  shown  in 
Figs.  19  and  20.  As  previously  noted,  these  two  variables 
with  high  correlation  coefficients  could  provide  the  neces¬ 
sary  link  tc  determine  the  mixed  layer  depth.  Since  DZ  and 
0Z BUZ  are  the  differences  in  depth  from  the  MLD  to  the  EOT 
and  EOZ,  respectively,  a  pair  of  simple  equations  define  the 
relationship . 
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MLD  =  EOT  ♦  DZ 

MLD  *  EOZ  +  EZBOZ 

note:  MID,  EOT  and  BUZ  are  negative  values  (depth) 

The  only  weak  "link”  in  the  two  equations  is  the  deter¬ 
mination  cf  the  EOT  and  BDZ  directly  from  SST.  To  improve 
the  correlation  coefficients  for  these  two  variables  as  sug¬ 
gested  by  the  regional  scatter  diagrams,  data  from  the  Sub¬ 
tropical  Hater  and  several  distant  data  points  which 
detracted  from  the  established  trends  were  removed.  In  sev¬ 
eral  cf  the  plots,  primarily  fcr  the  Temperature-dependent 
variables,  where  the  trend  was  actually  stronger  and  more 
accurate  with  data  from  the  Subtropical  Region  included,  the 
correlation  coefficients  were  improved.  The  correlation 
coefficients  fcr  the  best-fit  trends  are  listed  in  Table  v. 
Indications  are  made  to  show  which  coefficients  increased  or 
decreased  frcm  the  values  presented  in  Table  IV. 

An  example  cf  a  correlation  coefficient  that  decreased 
due  to  the  removal  of  ca*-.a  points  and  yet  provided  a 
stronger  trend  was  the  SLD  versus  DT5  combination  which 
changed  frcm  -0.484  to  -0.197.  In  this  case,  the  majority 
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Figure  20.  Regional  Scatter  Diagram  -  DZBUZ  Versu 


TABLE  V 

Best-Fit  Ccrcelati.cn  Coefficients  of  Thermal  Variables 


MU) 

.535y  SST 


-.197°  .27l'  DT5 
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of  data  points  were  grouped  about  the  DT5  value  of  0.4C  per 
5  d  and  a  HLE  value  of  -23  m.  The  attempt  to  capture  the 
trend  of  the  Transition  Region  was  made  by  excluding  the 
distant  scattered  points;  however,  this  detracted  from  the 
overall  trend  and  reduced  the  correlation  coefficient. 

A  large  number  of  correlation  coefficients  improved  sig¬ 
nificantly  by  this  method  of  data  exclusion.  Most  notice¬ 
able  is  the  variable  SST,  which,  (except  for  the  combination 
with  DT5)  increased  all  ether  coefficients  to  greater  than 
0.  645  with  six  increasing  above  0.900.  The  important  corre¬ 
lations  of  SST  with  BOT  and  BUZ  improved  from  -0.  203  and 
-0.682  to  0.776  and  -0.790,  respectively.  The  correlation 
coefficient  cf  SST  versus  EOT  even  reversed  sign. 

This  increased  strength  in  the  various  correlations 
between  variables  provided  the  confidence  to  infer  that  usa¬ 
ble  relationships  exist  in  the  vertical  thermal  profile. 
Since  the  variable  SST  had  the  greatest  number  of  strong 
correlations  with  the  other  variables,  the  statistical 
method  of  deriving  linear  regression  equations  using  only 
the  input  cf  SST  appeared  practical. 
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E.  HEGHESSICN  ECUAT ICN  ANALYSIS 

From  the  test  fit  correlation  coefficients  of  the  vari 
able  SST  in  combination  with  the  other  variables,  rhe  fcl 
lowing  linear  regression  equations  were  derived: 

MID  =  (1.2CC1)  SSI  -  43.300 

CIS  =  (0.03629)  SST  -  0.21573 

DT 15  =  (0.  16058)  SST  ♦  2.2271 

DT50  =  (0.48462)  SST  -  0.14105 

DT  10 0  =  (0.48506)  SST  ♦  0.  52999 

BCT  =  (3.9611)  SST  -  133.08 

TBOT  (0.61762)  SST  -  1.9317 

DT  =  (0.41729)  SST  ♦  0.94677 

DZ  =  (-2.7611)  SST  +  89.782 

E0Z  =  (-7.5247)  SST  -  0.00793 

TEUZ  =  (0.48510)  SST  -  0.66001 

DTBOZ  =  (0.55345)  SST  -  Q.U213Q 

DZEUZ  =  (7.7697)  SST  -27.820 


81 


To  construct  a  complete  vertical  temperature  prefile, 
since  the  variables  previously  mentioned  only  defined  the 
upper  domain  cf  the  thermal  structure,  a  new  variable  was 
needed  from  lower  depths.  The  temperature  at  460  m  was  cho¬ 
sen  because  it  was  the  lowest  datum  point  available  from  the 
XBT  data  provided  along  the  track.  This  type  of  variable 
was  not  particularly  desired  as  it  relates  a  temperature  to 
a  specific  depth  rather  than  to  a  feature  on  the  thermal 
profile.  A  better  variable  would  be  at  a  reference  depth 
where  the  temperature  gradient  becomes  isothermal  in  the 
deeper  regions  of  the  ocean.  However,  this  was  not  availa¬ 
ble  from  the  limited  depth  of  the  data  and  the  T460  variable 
was  established.  The  same  methods  of  analysis  were  employed 
on  this  variable  as  were  conducted  on  the  others.  When  com¬ 
pared  to  the  sea  surface  temperature,  a  correlation  coeffi¬ 
cient  of  0.831  was  computed  along  with  a  strong  linear 
trend  in  the  regional  scatter  diagram.  Therefore,  a  linear 
regression  equation  could  be  used  to  determine  the  bottom 
point  of  the  predicted  temperature  profiles  at  460  me-srs 
from  the  variable  SST.  The  linear  regression  equation  for 
T460  was  computed  as: 


T460  =  (0.21834)  SST  +  1.3109 


Two  ether  important  relationships  which  were  derived  by 
regression  analysis,  but  they  do  not  use  SST  as  the  indepen¬ 
dent  variable: 

DZ  =  (-0.85394)  BCT  -  12.693 
DZBUZ  =  (-1.0359)  BUZ  -  27.392 

These  equations  use  the  computed  values  of  30T  and  BUZ 
as  determined  from  the  SST  to  calculate  the  thickness  of  the 
seasonal  and  main  thermcclines .  The  magnitude  of  these 
thicknesses  (CZ,  DZEUZ)  when  added  to  the  BOT  or  BUZ, 
respectively,  can  then  define  the  mixed  layer  depth. 

Due  to  the  relatively  low  correlation  coefficient 
between  SST  and  MLD  for  the  SST-derived  MLD  equation 
(0.535),  possibly  a  combination  of  equations  which  define 
3LD  with  higher  correlations  could  improve  the  accuracy  in 
determining  the  mixed  layer  depth.  The  regression  equation 
for  the  SST-derived  BUZ,  with  a  correlation  coefficient  of 
-0.790  between  the  two  variables,  along  with  the  BUZ-aerived 
DZBUZ  mentioned  above  with  a  correlation  coefficient  of 
-0.983,  seemingly  provide  the  strongest  combination. 
Therefore,  the  group  of  equations  used  for  the  mixed  layer 
decth  determination  were: 


BUZ  =  (-7.5247)  S  ST  -  0.00793 

DZEUZ  =  (-1.0359)  BUZ  -  27.392 
MID  =  BUZ  ♦  DZ  BUZ 

Three  ctfcer  combinations  of  regression  equations  could 
have  been  used  to  determine  the  mixed  layer  depth,  but  were 
not  chosen  because  of  their  lower  correlation  coefficients 
between  variables.  The  equations  and  correlation  coeffi¬ 
cients  (r)  are  listed  belcw: 


BOT  *  (3.S611) 

SST  -  133.08 

r  =  0.776 

DZ  =  (-2.7611) 

S ST  ♦  89.732 

r  =  -0.645 

MLD  =  BCT  ♦  DZ 

BOT  =  (3.S6  11) 

SST  -  133.08 

r  =  0.776 

DZ  =  (-0.6  5  3  94)  EOT  -  12.693  r  =  -0.949 

MLD  *  BCT  ♦  DZ 
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EUZ 


(-7.5247)  SST  -  0.00793 


r  =  -0.790 


DZBOZ  =  (7.7697)  SST  -  27.820  r  =  0.766 
HLO  =  BtJZ  ♦  DZBUZ 

iith  the  MID  determined,  the  temperature  at  this  depth 
was  needed  tc  provide  the  fundamental  temperature  from  which 
the  thermccline  could  be  computed.  A  variable  TMLD  was 
added  to  represent  the  temperature  at  the  mixed  layer  depth 
and  was  defined  as  the  sum  of  TE0Z  and  DTBUZ.  Since  both 
temperature-dependent  variables  had  very  high  correlations 
with  the  SST,  0.949  and  0.960  respectively,  strong  confi¬ 
dence  could  be  placed  in  the  resulting  variable,  TMLD. 

An  additional  feature  in  calculating  TMLD  by  this  method 
was  to  allcw  for  a  slight  negative  temperature  gradient  tc 
occur  in  the  mixed  layer.  Many  ocean  thermal  structure  mod¬ 
els  assume  that  this  upper  layer  will  be  isothermal  when 
actually  a  realistic  profile  often  includes  a  temperature 
gradient . 

The  four  thermal  gradient  variables  DT5,  DT15,  DT50  and 
DT 100,  as  determined  from  the  SST,  were  subtracted  from  the 
temperature  and  depth  of  the  mixed  layer  to  form  the  thermo- 
cline.  This  structure,  along  with  the  SST  and  mixed  layer 
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depth  values,  provided  the  upper  portion  of  the  vertical 
temperature  prefile. 

In  the  region  below  the  seasonal  thermoclire  the 
temperature-dependent  variables,  1E0T  and  TBOZ,  and  depth- 
dependent  variables,  BOT  and  BOZ,  were  computed  directly 
from  the  input  cf  SST  and  applied  to  the  temperature  pro¬ 
file.  Several  times  during  the  plotting  of  these  points  in 
the  testirg  and  ccmpariscr  phase  of  the  thesis,  the  location 
of  the  DT100  variable  would  occur  below  both  the  BOT  and  BUZ 
or  between  them.  In  these  cases  it  was  simple  to  adjust  the 
sequerce  cf  pcints  to  allow  the  proper  construction  of  the 
temperature  profile. 

Finally,  in  the  generation  cf  the  temperature  profiles, 
the  variable  T46C  was  used  to  provide  the  lower  base  of  the 
structure.  Since  T460  was  an  artificially-generated  vari¬ 
able,  defined  at  a  constant  depth,  the  characteristic  small 
temperature  variations  at  that  depth  should  still  be  accu¬ 
rately  reflected. 

A  total  cf  nine  points  were  used  to  define  the  entire 
vertical  temperature  profile.  These  were  plotted  or.  a 
T2 KTKGNIX  618  graphic  display  system  and  connected  with  a 
rational  spline  method  that  allows  some  curvature  between 
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points.  This  was  done  to  portray  the  temperature  profiles 
of  the  region  acre  realistically. 

P.  THERMAL  EFCFILE  TEST  AND  COMPARISON  ANALYSIS 
1 .  TEST  1 

The  first  type  of  comparison  test  was  conducted  with 
tha  original  trade  of  data  used  in  the  generation  of  the 
predicted  thermal  profiles.  The  XET  temperature-depth  pro¬ 
files  of  the  Subarctic  Water  and  Transition  Zona  were  com¬ 
pared  with  the  computed  temperature  profiles  which  utilized 
the  sea  surface  temperature  of  the  actual  XET  trace  as  tha 
only  input.  Fig.  21  shows  the  comparison  between  the 
observed  and  predicted  thermal  structures  for  the  northern¬ 
most  XBT  station  on  the  track.  Appendix  F  lists  the  remain¬ 
ing  twenty  profile  comparisons  for  the  northern  regions 
excluding  the  Subtropic  Water. 

The  comparisons  generally  illustrated  a  strong  simi¬ 
larity  in  the  upper  layers.  The  mixed  layer  dapth  along 
with  the  slight  negative  temperature  gradient  in  the  layer 
were  realistically  predicted  for  a  majority  of  cases;  inter¬ 
nal  wave  metiens  are  assumed  responsible  for  the  largest 
variation  from  the  predicted  profiles.  The  upper  thermo- 
cline  region  closely  approximated  the  strong  thermal 


87 


•00  -490  -400  -990  -300  -390 


0  3-9  3  M  10  U  It  1M  30  22J 

TEMPERATURE  (DEGREES  O 


.  Observed  and  Predicted  Thermal  Profile 
Station  3  39. 


Figure  21 


gradient  measured  by  the  XBTs.  In  the  lower  thermcclins, 
greater  variability  occurred,  assumed  due  to  internal  waves 
and  thermal  advection,  resulting  in  poor  comparisons  as 
shown  in  Figs.  F.2,  F. 1 1  and  F. 12  for  XBT  stations  341,  350 
and  351. 

The  predicted  lower  regions  were  not  as  precise  when 
compared  to  the  observed  profiles  usually  due  tc  the  exis¬ 
tence  of  a  significant  temperature  inversion.  Lack  of  a 
model  for  this  feature,  which  was  not  attempted,  clearly 
affected  the  profiles  at  depths  greater  than  130  m.  Several 
predicted  profiles,  however,  did  attempt  to  reproduce  this 
structure  with  the  placement  of  the  DT100  variable  below  the 
BUZ.  This  was  shown  in  XET  stations  343-346  and  348  (Figs. 
F.  4  -  F.  7  and  F .  9) . 

Between  the  bottom  of  the  main  thermoclina  and  460 
m,  several  predicted  profiles  completely  failed  to  approxi¬ 
mate  the  observed  structure.  This  was  caused  by  the  lack  of 
variables  in  the  lower  layers  as  well  as  the  computer-gener¬ 
ated  curvature  between  data  points.  When  the  variables 
DT100  and  BUZ  were  very  close  together,  as  shown  for  XET 
stations  348  and  349  (Figs.  F.9  and  F.10)  ,  the  profile  would 
bulge  with  a  positive  or  negative  gradient. 
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Utilizing  only  nine  points  to  define  the  predicted 
thermal  profile  obviously  smooths  the  fine  structure  within 
the  theroocline  and  mixed  layer.  However,  due  to  the  small 
size  and  transient  nature  of  these  features,  their  removal 
would  not  severely  limit  the  usefulness  of  the  predicted 
structures . 

In  this  initial  test,  since  the  predicted 
temperature  structures  were  generated  and  compared  with  the 
same  XBT  data  from  the  meridional  track,  the  comparisons 
should  be  relatively  close.  The  predicted  profiles  would 
actually  represent  the  average  of  all  the  observed  thermal 
structures  as  the  variables  reflect  the  trends  along  the 
track. 

2.  1IST  2 

A  more  stringent  and  realistic  test  was  attempted 
that  defines  the  spatial  limitations  of  this  thermal  pre¬ 
diction  method.  Test  2  consists  of  comparing  observed  and 
predicted  temperature  profiles  which  occurred  during  the 
same  time  period  as  the  data  generation,  but  spatially 
removed.  A  possible  indication  as  to  a  limit  in  the  effec¬ 
tive  distance  would  be  the  observation  that  the  predicted 
prefile  begins  to  diverge  greatly  from  the  actual 


temperature  profile 


The  test  XBT  data  were  obtained  from  Fleet  Numerical 
Oceanography  Center  by  a  computer  search  for  all  available 
XBT  drops  which  occurred  in  September  from  1965  to  1980  in  a 
ren  degree  wide  strip  of  longitude  which  surrounds  the  orig¬ 
inal  track.  Inexplicably,  the  data  received  were  limited  to 
the  area  south  of  44N,  excluding  the  Subarctic  Water.  Lack 
of  time  prevented  an  attempt  to  acquire  these  data  whicn 
would  be  useful,  but  not  entirely  necessary,  for  analysis. 
The  temperature-depth  data  were  then  selected  by  noting  the 
positions  relative  to  the  original  data.  To  provide  an 
effective  test,  a  variety  of  locations  in  the  Transitional 
Zone  were  chosen.  However,  these  positions  were  constrained 
by  the  availability  of  data  in  the  region  and  appeared  as  a 
tonal  belt  perpendicular  tc  the  original  track.  Fig.  22 
shows  the  location  of  original  and  test  XBT  drops. 

The  temperature-depth  points  provided  by  FNOC  were 
then  plotted  and  connected  on  the  computer  graphical  display 
with  the  same  amount  of  curvature  applied  tc  the  profile  as 
before.  Ecth  observed  and  predicted  profiles  were  plotted 
together  again  for  analysis  and  comparison  purposes. 


Five  vertical  temperature  profiles  from  the  10-11 
September  period  were  compared  with  the  predicted  profiles 


as  generated  from  the  sea  surface  temperature  of  the 
observed  profiles.  The  test  produced  varied  results  with 
several  features  of  the  computed  temperature  profiles  com¬ 
paring  favorably  with  the  actual  thermal  structure.  The 
largest  differences  were  noted  for  the  two  profiles  located 
the  greatest  distance  from  the  meridional  track  at  145W. 
Figs.  23  and  24  illustrate  the  comparisons  at  these  loca¬ 
tions,  approximately  400  km  east  and  west  from  the  original 
track  of  data.  The  mixed  layer  depth  and  upper  thermocline 
regions  cf  the  predicted  profile  closely  resembled  the 
observed  profile,  but  both  structures  displayed  differences 
in  the  lower  regions. 

The  eastern  most  observed  profile  had  a  thin,  strong 
thermocline  which  gradually  weakened  with  depth.  This  con¬ 
trasted  markedly  with  the  strong  a*, d  sharply-defined  thermo¬ 
cline  of  the  predicted  profile.  The  regions  below  150  m 
displayed  similar  temperature  gradients  for  both  profiles. 

For  the  western  most  test  profile,  the  only  region 
to  differ  significantly  from  the  predicted  structure  was  the 
lower  domain  below  130  m.  The  location  of  the  two  close 
variables,  80Z  and  DT100,  along  with  the  computer-generated 
curvature  forced  the  predicted  profile  to  diverge  from  the 
weak  temperature  inversions. 
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The  three  other  test  cases,  shown  in  Figs.  25,  26 
and  27,  which  were  located  closer  to  the  original  track,  all 
displayed  strong  similarities  between  the  observed  and 
predicted  profiles.  Except  for  the  fins  structure  and  pos¬ 
sible  effects  of  internal  waves  cn  the  depth  of  the  thermo- 
clire,  the  computed  profiles  were  very  accurate.  The 
unusual  shape  of  the  observed  temperature  profile  at  the 
bottom  of  the  seasonal  theimocline,  in  Fig.  27,  and  close 
approximation  by  the  predicted  profile,  illustrates  the 
strong  ability  of  this  statistical  method  to  produce  thermal 
features  unique  to  a  certain  region. 

Tc  test  the  spatial  limitations  of  this  empirical 
method  further,  a  vertical  temperature  profile  taken  at 
Ocean  Weather  Station  "P"  during  the  data-ganeration  time 
period  was  examined.  Since  the  location  of  the  northern 
most  XBT  drop  cn  the  track  was  two  degrees  of  latitude  south 
of  OWS-P,  the  results  should  indicate  the  reliability  of  the 
methcd  outside  the  test  region.  Fig.  28  displays  the  com¬ 
parison  between  the  observed  and  predicted  profiles  for 
QWS-E. 

Frcm  the  surface  tc  the  middle  of  the  thermocline, 
the  two  profiles  were  very  similar,  with  the  mixed  layer 


96 


Figure  26.  Observe!  and  Predicted  Thermal  Profiles 
42. 3N,  144.  OW  For  102115Z  SSP  77. 


depth  accurately  depicted.  However,  the  lower  segment  of 
the  two  profiles  differed  by  2.0C.  With  the  actual  vertical 
temperature  prefile  being  considerably  cooler  than 
predicted,  a  different  regime  must  have  been  encountered. 
Tabata  (1961)  suggests  from  his  observations  in  the  region 
that  horizontal  transport  of  colder  water  from  the  vicinity 
of  the  Alaskan  Gyre,  lying  northwest  of  the  station,  may 
cause  a  decrease  in  temperature  during  the  summer  months. 

3.  TE§T  3 

The  third  test  included  both  spatial  and  time  dif¬ 
ferences  in  an  attempt  tc  forecast  the  vertical  thermal 
structure  up  to  three  weeks  in  the  future.  The  only  availa¬ 
ble  XBT  test  da*a  in  the  Transitional  Region  ware  located 
within  three  degrees  longitude  to  the  west  cf  the  original 
track  as  shewn  ir.  Fig.  29.  Vertical  temperature  profiles 
were  chosen  at  weekly  intervals  and  compared  with  the 
predicted  structures. 

At  the  end  cf  the  first  week,  the  observed  thermal 
profile,  as  shewn  in  Fig.  30,  displays  a  deeper  mixed  layer 
to  a  depth  of  4C  m  overlying  an  extremely  strong  negative 
temperature  gradient.  This  feature  was  not  represented  well 
by  the  predicted  profile  as  the  computed  MLD  and  thermocline 


variables  were  shallower  ir.  depth.  The  shape  of  the 
observed  nixed  layer  infers  strong  wind  and  convective  mix¬ 
ing,  so  surface  weather  charts  for  the  period  of  11-18  Sep¬ 
tember  1977  in  the  Northeast  Pacific  were  checked  fcr  the 
occurrence  of  a  significant  at  mo  spheric  event.  On  17  Sep¬ 
tember  1977,  a  cold  front  associated  with  a  low  pressure 
system  in  the  Gulf  of  Alaska  had  quickly  passed  through  the 
tes*  region  with  rain  and  35  knot  winds.  The  weather  during 
the  period  of  data  acquisition,  8-11  September  1S77,  as 
described  tsfcre,  was  mostly  cloudy  with  periods  of  light 
winds  and  scattered  showers  (FNOC,  1977).  These  contrasting 
atmospheric  situations  suggest  a  major  limitation  in  the 
statistical  model,  as  weather  effects  on  the  upper  vertical 
temperature  structure  must  be  properly  accounted  for  in  the 
predicted  profile. 

The  lower  half  of  the  thermccline  on  the  predicted 
temperature  prefile  closely  followed  the  observed  structure 
down  to  170  m.  Except  for  the  weak  temperature  inversion  at 
200  m,  which  was  not  included  ir.  the  empirical  model,  the 
lower  domain  behaved  well. 

The  next  test  profile  was  taken  or.  24  September 
1977,  thirteen  days  after  the  data  generation  for  the 
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predicted  structure.  Location  of  this  test  XBT  was  on  the 
southern  boundary  of  the  Transitional  Zone  approximately  200 
km  west  of  the  original  track.  The  predicted  temperature 
profile  in  Fig.  31  clearly  resembled  the  observed  vertical 
thermal  structure,  but  with  slightly  cooler  temperatures. 
Due  to  the  strong  winds  from  the  storm  a  week  earlier,  the 
observed  MID  remained  deep.  However,  a  transient  at  30 
depth  hac  formed  with  a  0.5C  increase  in  sea  surfac 
temperature  as  the  weather  in  the  region  during  the  secc.  _ 
week  was  dominated  by  a  strong  high  pressure  system.  The 
observed  profile*s  thermo cline  also  deepened  by  about  10  m, 
which  may  be  due  to  internal  waves,  or  it  may  have  been 
caused  by  warm  advection  from  the  south.  Immediately  below 
the  thermocline,  the  two  profiles  differed  slightly  because 
of  the  fine  structure  in  the  region.  The  predicted 
temperature  at  460  m  depth  remained  accurate  when  compared 
to  the  observed  temperature  throughout  ail  the  temporal 
tests. 

The  final  test  forecast  was  conducted  for  an  X3T 
location  209  km  to  the  west  of  the  north-south  track  on  30 
September  1977,  nearly  twc  weeks  after  the  basic  data  were 
acquired.  Another  cold  front  associated  with  a  low  pressure 
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system  moved  through  the  area  on  25-26  September  1977  and 
was  immediately  followed  t7  a  high  pressure  ridge.  The 
storm,  with  maximum  winds  of  40  knots,  maintained  the  40  m 
thick  isothermal  mixed  layer  as  shewn  in  Fig.  32.  The  wind 
mixing  had  increased  and  deepened  the  negative  temperature 
gradient  cf  the  thermocline  to  a  level  approximately  10  m 
below  the  predicted  profile.  In  the  lower  domain  of  the 
thermal  structure  a  large  difference  was  noted  due  to  the 
slight  temperature  inversion  of  the  observed  profile  and  to 
the  close  location  of  the  two  variables,  3UZ  and  DT100,  at 
130  m  on  the  predicted  prefile.  The  profiles  coincided  at 
the  460  m  depth  as  determined  by  the  variable  T460. 

The  accuracy  of  the  predicted  profile  after  three 
weeks  from  the  time  of  data  acquisition  for  the  statistical 
analysis  shows  an  overall  degradation.  However,  considering 
that  two  storm  systems  passed  through  the  region  and  the 
summer  heating  season  was  approaching  its  end,  the  fore¬ 
casted  profile  was  not  excessively  peer. 

G.  REVIEW  CF  PREDICTED  TEERMAL  PROFILE  RESULTS 

The  underlying  requirement  of  the  statistical  approach 
in  predicting  the  vertical  temperature  profile  is  to  acquire 
accurate  data  that  reflects  the  mean  thermal  conditions  for 
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a  particular 


s«ason  and  region.  Besides  the  acquisition 
techniques,  the  importance  in  precisely  defining  the 
regional  boundaries  is  essential.  Both  temperature  and 
salinity  should  be  used  to  establish  the  limits  of  usable 
data  when  applying  correlation  and  regression  analyses. 
Results  of  the  predicted  thermal  profiles  located  outside 
the  data  generation  area,  as  defined  by  a  certain  region, 
displayed  a  marked  degradation  in  the  lower  thermociine  and 
lower  domain  regions.  The  predicted  profiles  located  near 
the  area  generally  portrayed  the  observed  temperature  pro¬ 
files  realistically  within  the  spatial  limits  of  455  km 
about  the  original  data. 

The  results  cf  the  temporal  accuracy  of  the  predicted 
profiles  indicate  a  favorable  ability  to  forecast  the  ther¬ 
mal  structure  with  certain  limitations.  The  main  limiting 
factor  was  the  weather  which  must  be  included  in  the  analy¬ 
sis.  Clearly,  when  high  winds  pass  through  a  region  after 
the  vertical  temperature  data  is  acquired,  the  upper  thermal 
structure  will  be  altered.  Two  regions  of  the  thermal  pro¬ 
file  that  will  be  affected  are  the  depth  of  the  mixed  layer 
and  strength  of  the  thermociine.  For  short  periods  of  abcut 
one  week,  this  upper  domain  in  the  water  column  will 
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primarily  reflect  changes  due  to  the  wind  and  f: 


convec¬ 


tive  mixing.  *cr  time  scales  cf  greater  than  a  week,  an 
occasional  thermal  intrusion  from  the  Subarctic  or  Subtropic 
Region  may  be  an  additional  factor  that  could  alter  the 
thermal  structure.  Such  advective  processes  would  be  diffi¬ 
cult  to  forecast  by  the  statistical  method  in  regions  where 
large  meanders  and  eddies  occur,  but  they  should  net  greatly 
affect  a  relatively  quiescent  area  such  as  ir.  the  Tran¬ 
sitional  Zone. 

An  example  cf  another  problem  encountered  by  this  method 
of  thermal  structure  prediction  from  the  sea  surface 
temperature  is  the  process  of  heating  at  the  ocean's  sur¬ 
face,  called  the  "afternoon  effect."  Pig.  33  illustrates 
the  difference  in  the  predicted  profile's  structure  with  the 
observed  "afternoon  effect."  Since  the  sea  surface 
temperature  determines  the  thermal  profile  in  this  statisti¬ 
cal  model,  an  error  is  imposed  in  the  subsurface  structure, 
throughout  the  entire  depth.  With  the  surface  heating 
removed,  the  predicted  profile  would  closely  resemble  the 
observed  profile.  The  "afternoon  effect"  would  greatly 
hinder  the  ability  of  remote  sensing  to  determine  the  sub¬ 
surface  thermal  structure  by  this  method. 


1 10 


The  method  c£  determining  nixed  layer  depth  was  reviewed 
to  ensure  that  the  best  possible  value  was  computed  and  used 
ir.  the  empirical  scheme.  As  discussed  earlier,  the  ML D  pro¬ 
vides  a  base  for  the  variables  DT5,  DT15,  DT50,  DT100,  DT, 
DZ,  DTBUZ  and  DZBUZ  to  define  their  values  in  the  predicted 
temperature  profile.  The  ether  variables  BOT,  TBOT,  BUZ  and 
TBUZ,  besides  contributing  to  the  construction  of  the  pro¬ 
file,  provide  a  check  to  maintain  continuity  with  the  MLD- 
de terminal  variables.  The  technique  used  in  this  thesis  to 
compute  MLD  caused  no  problems  in  the  determination  of  the 
thermal  profile;  however  a  comparison  test  of  computing 
mixed  layer  depth  by  the  five  different  methods,  defined 
earlier  in  Chapter  IV,  was  conducted  to  ascertain  which 
method  is  the  mest  accurate. 

All  methods  were  provided  with  the  same  SST  values  from 
the  original  data  of  the  rcrth-scuth  track  and  compared  with 
the  actual  mixed  layer  depth.  The  computed  depths  were  then 
correlated  to  determine  the  best  technique.  Results  from 
the  test  indicated  all  methods  were  surprisingly  equal.  The 
methed  used  was  chosen  because  of  the  higher  correlations 
oetween  variables.  However,  the  simplist  method,  which  was 
the  SST-derived  nixed  layer  depth,  with  a  correlation  coef¬ 
ficient  of  0.535,  could  have  been  used  with  equal  success. 


V.  CONCLUSIONS  AND  RECOMMENDATIONS 

A.  CONCLUSIONS 

The  sea  surface  temperature  reflects  many  cf  rhe  physi¬ 
cal  processes  that  govern  the  oceans.  Use  of  r snore  sensing 
to  map  the  surface  thermal  structure  of  the  world’s  oceans 
accurately  and  rapidly  presently  exists,  but  the  determina¬ 
tion  of  precise  subsurface  vertical  temperature  structure 
from  these  surface  signatures  is  unrealized.  Possibly  only 
in  selected  ocean  regions  can  insight  be  gained  into  deter¬ 
mining  quantitative  relationships  for  directly  inferring  the 
subsurface  features  from  the  sea  surface  temperature.  In 
relatively  quiet  locations  with  uniformly  changing  proper¬ 
ties,  such  as  in  the  Transitional  Zone  of  the  Northeast 
Pacific,  many  strong  correlations  can  be  observed  within  the 
thermal  structure. 

The  results  cf  the  statistical  analysis  indicated  that 
the  strongest  relationships  within  the  vertical  temperature 
profile  occurred  with  the  thermally-dependent  variables 
(TBOT ,  Dl,  TEUZ,  DTBUZ)  as  they  varied  with  SST.  The 
depth-related  variables  (MID,  30T,  DZ,  BUZ,  DZ3UZ)  displayed 
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weaker  correlations  with  the  surface  temperature,  but  trends 
still  could  be  observed  through  the  thermal  structure.  The 
use  cf  scatter  diagrams  tc  assist  in  the  correlation  analy¬ 
sis  cf  variables  was  necessary  to  establish  regional  bounda¬ 
ries  and  the  test  fit  points  for  the  derivation  of 
regression  equations. 

An  important  region  within  the  vertical  temperature  pro¬ 
file  is  the  thermocline.  In  practical  terms,  its  depth, 
strength  and  gradient  strongly  influence  the  propagation  cf 
underwater  sound.  The  prediction  of  this  structure  by  use 
of  only  sea  surface  temperature  would  greatly  enhance  A  si 
capability.  Normally,  within  the  Northeast  Pacific,  the 
correlation  between  SST  and  thermocline  depth  is  not  very 
high;  however,  in  the  Transition  Zone  with  uniformly  chang¬ 
ing  structure  a  reasonable  predictability  can  be  obtained. 
The  thermocline  variables  (DT5,  DT15,  DT50,  DT103)  all  pro¬ 
vided  trends  from  which  SSI-related  values  were  computed  tc 
define  the  thermal  structure.  Overall,  the  most  precisely 
predicted  region  of  the  temperature  profile  was  the  thermo- 
cline.  The  depth  may  have  varied  due  to  internal  waves,  but 
the  negative  temperature  gradient  was  accurate  and  realistic 


when  compared  to  the  observed  structure. 


B.  RECOMMENEAUCNS 

In  future  research  cf  this  method,  an  expanded  sat  cf 
vertical  temperature  data  from  parallel  or  bisecting  tracks 
should  be  established  as  the  data  base  for  statistical  anal¬ 
ysis.  The  addition  of  an  east-west  track  of  data  would  pro¬ 
vide  latitudinal  trends  which  could  also  be  used  in  the 
predicted  temperature  profiles.  As  in  any  data  aguisiticn, 
a  consideration  as  to  the  minimum  amount  cf  data  points  *o 
determine  the  thermal  structure  statistically  should  be 
made.  Meanwhile,  a  number  of  no  less  than  20  profiles  ever 
a  distance  of  1000  km  with  spacing  between  points  not  to 
exceed  50  km  is  recommended. 

ether  suggestions  for  model  improvement  include  adding 
dynamical  and  thermodynamical  variables  such  as  wind  speed, 
direction  and  duration,  air  temperature,  radiation  flux,  and 
other  atmospheric  effects  which  could  alter  the  ocean  ther¬ 
mal  structure.  Most  of  these  factors  can  be  determined  by 
remote  sensing;  data  from  ether  satellite  sensors,  such  as 
an  altimeter,  SAP  and  Coastal  Zone  Color  Scanner  which  car- 
measure  sea  surface  signatures,  should  also  be  included. 

lo  determine  the  limitation  of  accuracy  for  the  statis¬ 
tically-determined  profile  method,  more  stringent  tests 


should  be  conducted  which  increase  the  spatial  and  rime  fac¬ 
tors.  In  certain  ocean  regions  one  could  possibly  maintain 
reasonable  accuracies  within  hundreds  of  kilometers  and  sev¬ 
eral  weeks  duration  from  the  acquisition  of  thermal  data. 
Different  ocean  locations  should  be  analyzed  and  nested  to 
determine  the  possible  regions  where  the  thermal  profiles 
vary  in  an  orderly  fashion  which  can  be  easily  modeled. 
Areas  where  mescscale  features,  such  as  eddies  and  fronts, 
are  constantly  altering  the  water  column  by  advective  pro¬ 
cesses  would  be  poor  candidates  for  this  method. 

Another  recommendation  is  to  test  in  other  seasons  of 
the  year  to  determine  if  the  statistical  approach  would 
apply.  The  difference  in  the  vertical  thermal  structure 
from  the  heating  to  cooling  season  is  sc  great  that  possibly 
there  could  be  only  a  certain  window  of  time  that  allows  the 
sea  surface  temperature  to  define  the  subsurface  structure. 

Finally,  a  method  fcr  determining  the  occurrence  of 
large  temperature  inversions  in  the  Subarctic  and  Transition 
Regions  is  needed.  This  important  thermal  feature,  because 
of  its  variation  in  magnitude,  depth  and  location,  affects 
greatly  the  propagation  of  sound  by  forming  a  sound  channel 
immediately  telcw  the  the rmoclin  e.  If  this  structure  can  be 
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BICMEEICAI  COMPOTES  PROGRAMS  (DIXON  AND  BROWN, 

1979) 

The  EICMEE  statistical  computer  programs  were  utilize! 
because  of  the  simple  and  flexible  programming  it  provided 
in  analyzing  the  thermal  data.  In  constructing  the  vertical 
temperature  profiles,  the  statistical  analysis  required  cor¬ 
relation  and  regression  equation  computation  and  scatter 
diagram  plotting  cf  the  thermal  variables. 

The  program  initially  selected  was  'Stepwise  Regression 
P 2R. '  This  program  computes  simple  statistics,  i.s., 
mean,  standard  deviation,  etc.,  correlation  coefficients  and 
linear  regression  equations  for  any  number  of  variables. 
The  statistics  used  in  the  computations  are  defined  as:. 


In  siiFl-  linear  regression  analysis,  where  there  is 
only  one  independent  variable,  the  coefficients  are  deter¬ 
mined  by  a  least  squares  method.  The  form  cf  the  equations 
are: 

y  *  A  ♦  Bx  (where  A  and  B  are  the  linear  coefficients) 

The  EICMEE  computer  program  used  to  print  the  regional 
scatter  diagrams  is  titled  ’Bivariate  (Scatter)  Plots 
P6 D. •  Regression  equations  for  regressing  variable  x  on  y 
and  y  on  x  can  be  requested  along  with  the  intersections  of 
the  regression  lines  with  the  axis  of  the  plot.  This  was 
useful  in  analyzing  the  data  to  determine  the  best  fit  for 
the  various  regions. 
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Figure  C.10.  Regional  Scatter  Diagram  -  SST  Versus  BOZ. 
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Figure  C.  13.  Pegional  Scatter  Diagram  -  53T  Versus  DZBUZ 
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Figure  D.a.  Fegicnal  Scatter  Diagram  -  SOT  Versus  DT100. 
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BEGICNAI  SCATTER  DIAGRAMS  -  BUZ  VERSUS  THERMAL 

VARIABLES 
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